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Abstract 


The  Affordable  Moving  Surface  Target  Engagement  (AMSTE)  project  attempts  to 
develop  affordable  solutions  to  the  precise  moving  target  surface  target  engagement 
problem.  Up  to  this  point,  most  of  the  error  analysis  performed  for  the  AMSTE  project 
has  been  at  the  error  variance  level,  generating  root-sum-square  (RSS)  total  errors  from 
error  budgets  consisting  of  constant  error  variances.  In  reality,  the  level  of  error  for  both 
Global  Positioning  System  (GPS)  positioning  and  radar  targeting  systems  is  highly 
dependent  upon  the  given  situation  (such  as  the  distance  between  sensor  and  target,  the 
altitude  differences,  etc.) 

This  research  generates  a  more  comprehensive  model  of  the  GPS  errors  based  upon 
the  underlying  physics  of  the  situation.  It  focuses  on  differential  tropospheric  errors  and 
multipath,  as  these  are  the  primary  error  source  in  a  differential  GPS  targeting  system. 

In  addition  to  the  error  model  development,  a  code-based  differential  GPS  and 
differential  ranging  approach  is  implemented  in  simulation  using  a  Kalman  filter.  This 
approach  uses  GPS  measurements  collected  by  each  of  the  sensors  and  the  weapon,  and  it 
uses  ranging  measurements  from  the  sensors  to  the  bomb  and  the  target.  Multiple  cases 
are  run  varying  1)  the  number  of  GPS  satellite  measurements  tracked  by  each  receiver,  2) 
whether  or  not  the  common  GPS  errors  are  estimated,  and  3)  whether  or  not  the  bomb  is 
tracked  with  the  same  radar  sensors  that  are  tracking  the  target.  The  horizontal  DRMS 
position  error  during  the  terminal  phase  of  the  bomb  trajectory  drops  from  about  6  meters 


to  about  3.5  meters. 


USING  GPS  AS  A  REFERENCE  SYSTEM  TO  HIT  A  MOVING  TARGET 


1.  Introduction 

Over  time,  the  method  of  damaging  an  enemy’s  surface  infrastructure  or  hardware 
has  changed  dramatically.  In  the  earliest  days  of  cannons,  two  or  three  test  shots  were 
required  to  determine  the  range  to  the  targets.  Environmental  and  technical  factors  such 
as  weather  and  poor  targeting  repeatability  made  the  firing  of  cannons  a  guessing  game. 

In  today’s  warfighter  environment,  the  ability  to  engage  and  destroy  fixed  or  stationary 
surface  targets  has  been  honed  due  to  modem  technology.  As  adversaries  begin  to  adopt 
mobility  as  a  means  of  survival,  an  affordable  and  precise  moving  surface  target 
engagement  capability  will  be  necessary.  Currently  the  use  of  precise  guided  munitions 
(PGMs)  using  laser-guidance  with  a  “man-in-the-loop”  has  some  benefits  for  mobile 
targeting  [44].  In  the  effort  to  redefine  close  air  support  (CAS),  there  is  the  desire  to  have 
the  ability  to  transmit  updated  information  about  moving  targets  to  a  weapon  in  flight 
[41].  The  Affordable  Moving  Surface  Target  Engagement  (AMSTE)  program’s  goal  is 
to  provide  such  a  capability. 

In  the  development  of  any  network  or  weapon  system,  there  is  a  need  to  model 
expected  errors.  To  this  point,  most  of  the  error  analysis  performed  for  the  AMSTE 
program  has  been  at  the  error  variance  level,  generating  root-sum-square  (RSS)  total 
errors  from  error  budgets  consisting  of  constant  error  variances.  In  reality,  the  level  of 
error  for  GPS  positioning  and  targeting  systems  is  highly  dependent  upon  the  current 
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environment  in  which  it  is  employed  (i.e.,  distance  between  sensor  and  target,  altitude, 
time  of  day,  etc.) 

The  purpose  of  this  study  is  to  generate  more  comprehensive  models,  primarily  of 
the  GPS  errors,  and  look  at  the  sensor  errors,  based  upon  the  underlying  physics  of  the 
situation.  Differential  tropospheric  and  multipath  errors  will  be  the  focus,  since  these  are 
the  primary  error  sources  for  a  differential  GPS  targeting  system  to  be  employed  by  the 
AMSTE  program. 

1.1.  Background 

In  this  section,  the  Global  Positioning  System  (GPS)  will  be  briefly  covered. 
Additionally,  the  purpose  and  goals  of  the  AMSTE  program  will  be  discussed. 

1.1.1.  Global  Positioning  System  Overview 

The  Global  Positioning  System  (GPS)  is  a  satellite-based  radio-navigation  system 
established  by  the  U.S.  Department  of  Defense  for  military  positioning  applications  and, 
as  a  by-product,  it  is  also  serving  the  civilian  community.  The  system  provides  accurate, 
continuous,  worldwide,  three-dimensional  position,  velocity,  and  time  information  to 
users  with  the  appropriate  receiving  equipment. 

GPS  is  broken  into  three  segments:  the  space  segment,  the  control  segment  and 
the  user  segment.  The  space  segment  is  the  constellation  nominally  consisting  of  24 
satellites  arranged  in  6  orbital  planes  with  4  satellites  per  plane.  Each  satellites  orbit  has 
a  period  of  approximately  12  hours  and  a  55°  inclination  angle  with  respect  to  the 
equatorial  plane.  The  control  segment  is  a  worldwide  ground  control/monitoring  network 
monitoring  the  health  and  status  of  the  satellites.  These  monitor  stations  measure  signals 
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from  the  satellites  (SVs)  that  are  incorporated  into  orbital  models  for  each  satellite.  The 
models  compute  precise  orbital  data  (ephemeris)  and  SV  clock  corrections  for  each 
satellite.  The  Master  Control  Station  uploads  ephemeris  and  clock  data  to  the  SVs.  The 
Master  Control  facility  is  located  at  Schriever  Air  Force  Base  in  Colorado.  The  third 
segment  is  the  user  segment.  This  consists  of  all  GPS  receivers.  The  receivers  convert 
the  satellite  signals  into  position,  velocity,  and  time  estimates  [19]. 

Currently,  the  satellites  transmit  two  carrier  signals.  The  LI  frequency  (1575.42 
MHz)  carries  the  navigation  message.  Coarse  Acquisition  code  (C/A  code),  and  Precise 
code  (P-Code)  signals.  This  frequency  is  available  to  civil  users  as  the  Standard 
Positioning  Service  (SPS).  The  L2  frequency  is  modulated  by  the  P-Code.  Authorized 
users  with  cryptographic  equipment  and  keys  along  with  specially  equipped  receivers  use 
the  L2  frequency,  part  of  the  Precise  Positioning  Service  (PPS). 

The  navigation  message  is  a  50  Hz  signal  consisting  of  data  bits  that  describe 
parameters  such  as  the  GPS  week,  range  accuracy  prediction,  satellite  health,  clock 
corrections,  and  broadcast  ephemeredes.  The  C/A  code  is  a  repeating  1.023  MHz  Pseudo 
Random  Noise  (PRN)  Code.  This  “noise-like”  code  spreads  the  spectrum  over 
approximately  a  IMHz  bandwidth.  It  repeats  every  1023  bits  (one  millisecond).  Each 
satellite  has  its  own  unique  PRN  code.  The  C/A  code  on  the  LI  carrier  is  the  basis  for  the 
civil  SPS.  The  P-Code  is  a  very  long  (seven  days)  10.23  MHz  PRN  code.  In  the  Anti- 
Spoofing  (AS)  mode  of  operation,  the  P-Code  is  encrypted  into  the  Y-Code.  The 
encrypted  Y-Code  requires  a  classified  AS  Module  for  each  receiver  channel.  This  is  for 
use  only  by  authorized  users  with  cryptographic  keys.  This  P  (Y)-Code  is  the  basis  for 
the  PPS.  Navigating  by  the  pseudo-ranges  can  be  described  in  Figurel. 
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The  GPS  Navigation  Solution 

The  estimated  ranges  to  each  satellite  infersect  within  a  small  region  when  the  receiver  clock  bias  is 
comedy  estimaied  and  added  to  each  meacuied  relative  range. 

- -  MToana  oas 

Figure  1.  Illustrated  Description  of  GPS  Navigation  Solution  [9] 

The  position  of  the  receiver  is  where  the  pseudo-ranges  from  a  set  of  satellites 
intersect.  It  is  determined  from  multiple  pseudo-range  measurements  at  a  single 
measurement  epoch.  The  pseudorange  measurements  are  used  together  with  satellite 
position  estimates  based  on  the  ephemeris  data  sent  by  the  satellites.  This  orbital  data 
allows  the  receiver  to  compute  the  satellites’  positions  in  three  dimensions  the  instant 
they  sent  their  respective  signals.  A  minimum  of  four  satellites  is  required  to  determine 
the  three  position  dimensions  and  time.  Time  is  used  to  correct  the  offset  in  the  receiver 
clock,  allowing  receiver  clocks  to  be  less  expensive  and  less  accurate.  Using  five  or 
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more  satellites  provides  redundancy,  greater  fix  certainty,  and  can  allow  detection  of  out- 
of-tolerance  signals  [9]. 

1.1.2.  AMSTE  Program 

Existing  technology  (e.g.,  sensors,  communications,  and  weaponry)  supports 
moving  surface  target  engagement,  but  generally  requires  sophisticated  seeker  systems, 
humans-in-the-loop,  or  dispersive  area-effect  munitions  [43].  These  approaches  can 
result  in  expensive  weapon  systems,  high  risk  to  human  life,  and  the  potential  for 
collateral  damage  to  unintended  vehicles[43].  The  goal  of  the  AMSTE  program  is  to 
develop  and  demonstrate  a  new  strike  capability:  the  ability  to  target  moving  surface 
threats  from  long  range  and  to  engage  rapidly  those  threats  with  precision,  stand-off 
weapons.  This  capability  will  ultimately  enable  a  robust,  dynamically-controlled  system 
to  engage  multiple  dispersed  mobile  targets  accurately  in  complex  traffic  and  all  weather 
conditions  [6].  The  key  areas  are  described  as: 

•  Dynamic  -  contains  the  inherent  flexibility  to  create  a  custom  solution  for 
each  engagement 

•  Robust  -  degrades  gracefully  with  dynamic  changes 

•  Accurate  -  provides  precision  tracking  and  weapon  guidance  to  destroy 
the  moving  or  stationary  target  of  interest 

•  Multiple  -  has  the  ability  to  plan,  maintain  track,  and  conduct  parallel 
engagements 

•  Dispersed  -  has  flexibility  to  operate  over  an  entire  theater 

•  Mobile  targets  -  robust  to  target  move-stop-move  behavior  cycles 

•  Complex  traffic  -  maintains  target  track  in  dense  background  traffic 
environment 
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•  All  weather  -  provides  a  targeting  solution  in  any  weather  condition 

The  AMSTE  program  investigates  the  concept  of  leveraging  recent  advances  in 
sensor  technology  to  provide  an  affordable  solution  to  precise  moving  surface  target 
engagement  [1].  The  fundamental  concept  investigated  is  to  use  a  network  of  ground 
moving  target  indicating  (GMTI)  radar  and  synthetic  aperture  radar  (SAR)  systems  to 
provide  a  precision  fire  control  tracking  solution  on  moving  and  intermittent  surface 
targets.  These  sensors  provide  all-weather  capability,  and  the  use  of  multiple  GMTI 
sensors  provides  the  required  precision  targeting  and  will  be  provided  to  update  low-cost, 
precision  guided  munitions  in-flight  for  precise  engagement  of  a  moving  surface  target 
[2]. 

Current  studies  have  shown  that  obtaining  the  required  level  of  accuracy  is 
possible  with  expected  advancements  [11].  The  real  challenge  is  maintaining  a  precise 
estimate  of  the  location  of  the  target  [1]. 

1.2.  Problem  Statement 

The  problem  for  this  research  are  to  1)  identify  GPS  errors  that  are  significant  for 
the  AMSTE  program,  and  2)  determine,  through  simulation,  the  effects  of  GPS  and  radar 
errors  on  the  navigation  solution,  for  a  variety  of  implementation  schemes.  Additionally, 
the  simulation  results  will  be  used  to  demonstrate  the  benefit  of  tracking  both  the 
munition  and  target  with  the  same  sensor. 

1.3.  Scope 

Although  the  AMSTE  program  has  a  vested  interest  in  the  development  and 
accuracies  of  the  GMTFSAR  sensors,  this  thesis  focuses  on  GPS  and  its  contributing 
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errors.  The  GMTI/SAR  sensors  will  be  modeled  using  a  fairly  straightforward 
range/range-rate  representation,  as  described  in  Chapter  3. 

1.3.1.  Assumptions 

Typical  of  any  simulation,  assumptions  are  made  in  this  thesis  to  reduce  the 
complexity  of  the  development,  design,  and  analysis  of  the  GPS-based  system 
environment  model  and  Kalman  filter  simulations. 

1.  There  is  only  one  target  that  will  be  in  view  of  the  system  sensors  throughout  the 
scenario. 

2.  While  start  and  stopping  and  hiding  in  the  terrain  (e.g.  forests,  mountains)  would 
be  an  excellent  ploy  to  counteract  a  moving  target  system,  this  thesis  limits  the 
target  to  an  environment  similar  to  a  desert  where  the  vehicle  will  always  be 
moving. 

3.  The  location  of  the  Kalman  filter  processing  is  not  a  concern,  and  the  data  from 
all  sensors  is  sent  without  computation  or  transmission  delays. 

4.  Minimizing  the  error  in  the  relative  distance  between  munition  and  target  is  the 
focus.  The  absolute  position  of  the  munition  and  the  target  is  not  critical. 

5.  The  receivers  on  the  platforms  and  munition  are  all  identical.  While  each 
receiver  will  have  different  multipath  and  receiver  noise,  they  are  the  same  type, 
so  their  clocks  will  be  modeled  with  the  same  error  parameters. 

6.  Only  GPS  pseudorange  measurements  will  be  used.  While  using  carrier 
smoothed  or  phase  measurements  would  lead' to  more  accurate  results,  these 
procedures  are  outside  the  scope  of  this  thesis. 
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7.  Dual  frequency  receivers  will  be  used.  This  neglects  the  ionosphere  error  which 
will  be  discussed  in  Chapter  2. 

8.  Except  for  measurement  noise,  the  GMTI/SAR  sensors  are  identical. 

1.4.  Methodology 

The  premise  of  this  thesis  lies  in  the  modeling  and  building  of  a  true  environment 
model  and  a  system  upon  which  a  Kalman  filter  is  to  be  based.  The  environment  will  be 
built  using  the  true  ranges  and  rates  of  all  the  satellites,  platforms,  munition  and  target. 
Appropriate  errors  will  be  added  to  the  true  ranges  and  rates  to  generate  simulated 
measurements.  These  simulated  measurements  will  be  used  by  the  algorithm,  which  is 
essentially  a  Kalman  filter,  to  estimate  the  position  and  velocities  of  the  sensors,  bomb. 


and  target. 


2.  Background 


2.1.  Overview 

This  chapter  provides  background  that  is  needed  to  understand  the  true  environment 
model  and  system  Kalman  filter.  It  will  describe  differential  GPS,  the  errors  associated 
with  the  GPS  system,  how  differential  GPS  affects  those  errors,  and  errors  related  to  SAR 
range  and  range  rate.  The  chapter  closes  by  describing  the  basis  of  the  extended  Kalman 
Filter. 


2.2.  Differential  GPS  (DGPS) 

The  differential  GPS  concept  is  used  to  enhance  standalone  GPS  accuracy  by 
removing  common  or  correlated  errors  from  two  or  more  receivers  viewing  the  same 


satellites.  Figure  2  gives  a  graphical  representation  of  DGPS. 


Satellite 

“a” 


(Reference 

Receiver) 


Figure  2.  Differential  GPS  Concept 
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In  the  basic  form  of  DGPS,  one  of  these  receivers  is  called  the  reference  receiver. 
This  receiver’s  precise  position  is  known.  Because  of  this,  it  can  determine  the  errors  in 
the  pseudorange  measurement  from  each  satellite  in  view  by  subtracting  the  pseudorange 
from  the  satellite-to-reference  station  geometric  range.  Some  of  these  errors  would  be 
common  to  other  nearby  receivers  (atmospheric  delays,  satellite  clock).  The  reference 
station  would  then  send  out  those  errors,  called  differential  corrections,  to  the  mobile 
receiver  which  subtracts  these  common  errors  from  their  pseudoranges.  This  results  in  a 
more  accurate  position  solution,  on  the  order  of  l-4m  (Icr)  with  a  code-only  solution  [19, 
30, 43]. 

There  are  many  different  ways  to  implement  DGPS,  leading  to  different  levels  of 
accuracy.  As  mentioned  above,  using  code-only  measurements  (the  easiest  to 
implement)  has  an  accuracy  of  l-4m  (la).  Carrier-smoothed  code,  which  combines  the 
carrier-phase  precision  with  the  code  absolute  (but  noisy)  measurement,  can  have  an 
accuracy  of  0.1-0.5  m  (la)  [30, 34].  A  third  technique  is  to  use  the  carrier-phase 
measurement  and  resolve  the  carrier-phase  integer  ambiguities.  While  carrier-phase 
ambiguity  resolution  yields  an  accuracy  on  the  centimeter  level  (la),  carrier-phase  is  the 
most  difficult  to  implement  because  of  the  need  to  determine  the  integer  ambiguity  in  the 
measurement  [19,  21,  30]. 

For  this  study,  only  code  (i.e.,  pseudorange)  measurements  are  used.  There  are  a  few 
benefits  to  using  code  measurements:  1)  military  receivers  generally  only  give 
pseudorange  measurements,  and  2)  it  is  the  best  place  to  start  because  it  is  the  simplest 
form  and  easy  to  implement. 
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2.3. 


GPS  Errors 


Figure  3.  GPS  Errors  [12] 


Figure  3  shows  some  of  the  GPS  errors  that  affect  the  signal  from  the  satellite  to  the 
receiver.  GPS  errors  are  usually  classified  into  three  areas:  satellite,  atmospheric,  and 
receiver  errors.  These  errors  along  with  the  true  range  between  the  satellite  and  receiver 
make  up  the  pseudorange,  as  shown  in  the  following  equation: 

yo  =  r  +  c(<J  t„  -  t,,  +  t^„p  +  S  ti„„„  +  S  ^  t„,p)  (2-1) 

where: 


r 

=  geometric  range 

c 

=  speed  of  light 

=  receiver  clock  error 

=  satellite  clock  error 

^teph 

=  error  due  to  broadcast  ephemeris 

=  delay  due  to  the  troposphere 

=  delay  due  to  the  ionosphere 

^  ^  noise 

=  receiver  noise 

=  multipath  error 

Each  of  these  errors  will  now  be  discussed  in  detail. 
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2.4.  Satellite  Errors 


There  are  two  errors  associated  with  the  satellite.  These  are  satellite  clock  and 
satellite  ephemeris  errors. 

2.4.1.  Satellite  Clock 

GPS  satellites  use  atomic  clocks  to  control  the  onboard  timing  operations,  including 
the  generation  of  the  broadcast  signal.  Atomic  clocks,  which  are  typically  cesium  or 
rubidium,  are  inherently  stable.  Nevertheless,  they  do  drift,  and  over  time  the  5t  can  be  as 
big  as  1msec  [19].  This  can  translate  up  to  a  300-km  pseudorange  error.  Cloek 
correction  parameters  are  sent  in  the  navigation  message  and  are  implemented  by  the 
receiver  to  produce  estimates  of  the  actual  clock  errors.  The  correction  parameters  are 
implemented  by  the  receiver  using  the  following  second-order  polynomial: 

(S'  —  Ufo  +  afj(t  -  tgc)  ^f2(t  ~  toe)  (2"2) 

where 

Ufo  =  clock  bias 

Ufi  =  clock  drift 

fl/2  =  frequency  drift 

toe  =  clock  data  reference  time 

t  =  current  time  epoch 

Atr  =  correction  due  to  relativistic  effects 

Although  the  clock  error  estimate  removes  most  of  the  error,  some  residual  error 
remains,  and  this  residual  error  is  what  contributes  to  ranging  errors  in  the  pseudorange 
measurements.  The  resulting  ranging  errors  induced  by  the  satellite  clock  are  typieally 
on  the  order  of  3.0m  (la) [19]. 
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2.4.2.  Ephemeris  Errors  (  S 

Ephemeris  errors  result  when  the  GPS  message  does  not  transmit  the  correct  satellite 
location.  Estimates  of  ephemeredes,  as  in  the  clock  corrections,  are  uplinked  to  the 
satellite  and  broadcast  in  the  navigation  message.  Typically,  the  radial  component  of  this 
error  is  the  smallest.  The  tangential  and  cross-track  errors  may  be  larger  by  an  order  of 
magnitude.  Fortunately,  the  larger  components  do  not  affect  ranging  accuracy  to  the  same 
degree  as  the  radial  component  [36].  Generally,  the  ephemeris  errors  result  in  a  3m  (Ict) 
in  the  pseudorange  [30].  However  if  post-processing  is  an  option,  precise  ephemeris  data 
can  be  used.  Precise  orbits  of  the  satellites  are  calculated  using  days  of  data  from 
hundreds  of  reference  stations.  Precise  orbits  are  available  from  various  sources  (such  as 
National  Geodetic  Survey),  and  they  are  typically  accurate  to  within  6  cm  (la)  [30]. 

2.5.  Atmospheric  Errors 

If  a  vacuum  existed  between  the  satellites  and  the  receiver,  there  would  not  be  any 
error  to  due  the  propagation  of  the  signal.  The  atmosphere,  however,  affects  the 
transmission  and  the  reception  of  the  signal.  The  effects  of  the  ionosphere  and  the 
troposphere  are  normally  treated  as  separate  errors. 

2.5.1.  Ionosphere  (S 

The  ionosphere  is  the  part  of  the  atmosphere  consisting  of  free  electrons  that  ranges 
from  approximately  50km  to  1000km  above  the  surface  of  the  Earth  [19,22,  36].  It  is  a 
result  of  the  ultraviolet  radiation  from  the  sun  splitting  the  atmospheric  molecules  into 
ions  and  free  electrons.  The  density  of  these  electrons  determines  the  magnitude  of  effect 
on  the  GPS  signal.  Electron  density  is  normally  quantified  by  counting  the  number  of 
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electrons  in  a  vertical  column  with  a  cross-sectional  area  of  Im^.  This  is  called  the  Total 
Electron  Content  (TEC).  Figure  4  gives  an  illustration  of  the  TEC  [34]. 


Ionospheric: 

shell 


earth 


1  meter  cross- 
sectional 
area  tube 

. .v.,..:  • 

■  ■■ 


one  TEC  =  10^^  electrons  in 
column 

one  TEC  =  0.5  meters  delay 


Figure  4.  Ionosphere  Total  Electron  Content  (TEC)  Unit. 


TEC  (units  of  electrons/m^)  varies  with  time  of  day,  user  location,  satellite  elevation 
angle,  season,  ionizing  flux,  magnetic  activity,  sunspot  cycle,  and  scintillation  [19]. 
The  amount  of  delay  of  the  group  (pseudorange)  can  be  expressed  as: 


iono,g 


40.3  TEC 


(2-3) 


where 


^Siono,g  =  group  (pseudorange)  delay  (m) 

/  =  carrier  frequency  (LI  or  L2  for  GPS) 

The  delay  of  the  carrier-phase  delay  is 

40.3  TEC 


(2-4) 


The  ionosphere  advances  the  carrier-phase  by  the  same  magnitude  of  the  group  delay. 
Due  to  this,  a  dual  frequency  receiver  can  estimate  the  ionospheric  delay  and  remove  it. 
The  LI  ionospheric  delay  is  calculated  by 
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=  LI  ionospheric  delay  (m) 

/, ,  /^  =  LI  and  L2  carrier  frequencies 

Pn  ’  Phi  -  LI  and  L2  pseudorange  measurements 

The  L2  delay  can  be  calculated  by 


A5 


iono,corri^2 


2 


AS 


iono,corri^ 


(2-5) 


(2-6) 


While  applying  these  AS  delay  terms  completely  removes  the  ionospheric  error, 
multipath  and  noise  are  still  present,  adding  some  additional  error  [19,30].  For  this 
thesis,  it  is  assumed  that  dual  frequency  receivers  have  removed  the  ionospheric  delay. 
2.5.2.  Troposphere  ( S  ) 

The  troposphere  is  the  atmosphere  that  extends  from  the  ground  to  approximately  10 
km  above[19, 47,48].  Unlike  ionosphere,  the  tropospheric  delay  is  not  frequency 
dependent.  Within  the  range  of  the  GPS  frequencies,  the  troposphere  delays  both  the 
code  and  carrier  observations.  This  delay  is  dependent  on  the  tropospheric  refractive 
index,  which  is  a  function  of  the  local  temperature,  pressure,  and  relative  humidity. 

The  tropospheric  effect  can  be  divided  into  a  hydrostatic  (dry)  delay  and  a  wet  delay. 
The  hydrostatic  delay,  which  is  the  larger  effect  (90  to  95%  of  the  total  refraction),  is 
caused  primarily  by  N2  and  O2  molecules  [36].  The  hydrostatic  delay  at  zenith  (i.e., 
straight  up)  is  normally  around  2.3m,  and  it  varies  with  local  temperature  and 
atmospheric  pressure  in  a  reasonably  predictable  manner  [36, 47] .  Due  to  its 


predictability,  it  can  be  modeled  and  removed  with  an  accuracy  of  a  few  millimeters  or 
better  using  a  surface  model  (including  pressure,  temperature,  and  humidity)  [47]. 

The  second  part,  the  wet  delay,  is  generally  smaller,  with  delays  of  1-80  cm  at  zenith 
(depending  on  humidity  levels)  [36].  However,  since  the  wet  component  can  vary  10- 
20%  in  a  few  hours  [36],  it  is  more  difficult  to  model  and  remove  based  on  standard 
tropospheric  models  using  surface  measurements  [36,47].  The  residual  tropospheric 
delay  remaining  after  the  applied  model  is  mostly  due  to  the  wet  component. 

The  troposphere  delay  can  be  expressed  as 

^tropo=\i^-^)ds  (2-7) 


or 


AS,..,..  =10-*/^* 


(2-8) 


where 

n  =  index  of  refraction 

N  =  refractivity 

The  refractivity  N  can  be  divided  into  hydrostatic  and  wet  components.  Hence  the  above 
equation  can  be  written  as 

=  lO^J  N„,,,,ds  +  \0^\ N„„ds  (2.9) 


or  symbolically. 


d _ —  d  .  +  d 


tropo 


Hydro 


^Wet 


(2-10) 


where 

^ tropo  -  tropospheric  delay 
d Hydro  =  hydrostatic  delay 
‘^wei  =  wet  delay 
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Propagation  delays  at  various  elevation  angles  are  determined  from  the  zenith  delays 
and  mapping  functions.  As  the  zenith  delay  can  be  expressed  as  the  sum  of  the 
hydrostatic  and  wet  components,  mapping  functions  can  be  developed  in  order  to  map 
separately  the  hydrostatic  and  wet  components  [48].  In  general  the  tropospheric  delay  is 
represented  as 

dtrcpo  =  X  “Hy™,  (e)  +  X  ffV„  (f )  (2-11) 

where 

^ Hydro  =  hydrostatic  zenith  delay 

=  wet  zenith  delay 
=  hydrostatic  mapping  function 
^wet  (^)  “  mapping  function 
e  =  elevation  angle 

Much  research  has  gone  into  the  creation  and  testing  of  tropospheric  models  to 
compute  the  refractivity  N  along  the  path  of  signal  travel  [8,10,27,33,35,47,48].  The 
various  tropospheric  models  differ  primarily  with  respect  to  the  assumptions  made 
regarding  the  vertical  refractivity  profiles  and  the  mapping  of  the  vertical  delay  with 
elevation  angles  [56].  Even  though  a  tropospheric  model  removes  most  of  the  error,  there 
is  almost  always  a  residual  unmodeled  error. 

While  there  are  many  suitable  models  available,  this  research,  employed  the  modified 
Hopfield  model  because  it  has  proven  to  be  a  reliable  model  for  a  wide  variety  of 
environments.  The  thorough  explanation  of  this  model  is  available  in  [56]. 
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2.6.  Receiver  Errors 


While  there  are  various  types  and  models  of  receivers,  all  are  affected  by  receiver 
clock  errors,  multipath  errors,  and  receiver  noise.  Each  of  these  receiver-based  errors 
will  be  described  in  the  sections  that  follow. 

2.6.1.  Receiver  Clock  (St^  ) 

Similar  to  satellite  clock  errors,  receiver  clock  error  is  caused  by  the  inaccuracies  in 
the  receiver  clock.  However  unlike  the  satellite  clocks,  which  use  atomic  oscillators,  the 
receivers  typically  use  quartz  crystals,  which  have  much  worse  long-term  frequency 
stability.  Due  to  this  poor  stability,  receiver  clock  errors  are  typically  large,  but  this  error 
is  explicitly  estimated  as  part  of  the  navigation  solution. 

2.6.2.  Multipath  ( 6 

Multipath,  the  phenomena  whereby  a  signal  travels  from  a  transmitter  to  a  receiver 
via  multiple  paths,  is  due  to  reflection  and  diffraction.  It  is  a  major  source  of  error 
remaining  when  using  differential  GPS,  because  it  is  uncorrelated  between  receivers,  so  it 
does  not  cancel  out  [4,15,  28,32,35].  It  distorts  the  signal  modulation  (code)  and 
degrades  accuracy  in  conventional  and  differential  systems.  Multipath  also  distorts  the 
phase  of  the  carrier  and  hence  degrades  the  accuracy  of  the  carrier-phase-based  systems 
[3,4].  Additionally,  because  carrier-phase  algorithms  often  employ  pseudorange 
measurements  for  initialization  (ambiguity  resolution)  purposes,  multipath  contamination 
of  the  pseudorange  can  increase  the  time  required  for  initialization  [3,4].  Multipath- 
induced  errors  are  more  difficult  to  quantify,  because  they  are  a  function  of  the 
environment  in  which  the  particular  GPS  receiver  is  operating,  as  well  as  the  receiver 
design  [3,4]. 
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Multipath  has  been  the  subject  of  much  research  and  can  be  characterized  by  four 
basic  parameters: 

S  :  the  relative  time  delay  between  the  direct  and  reflected  signals 
a  :  the  amplitude  of  the  multipath  signal  relative  to  that  of  the  direct  signal 
(p  :  the  phase  of  the  multipath  signal  relative  to  that  of  the  direct  signal 
0  :  the  rate  of  change  of  the  phase  of  the  multipath  signal  relative  to  that  of  the 
direct  signal 


Figure  5  gives  an  illustration  of  the  relationship  between  the  relative  phase  and  delay  of 


multipath  and  the  resultant  range  error. 


Figure  5.  Multipath-Induced  Error  as  a  Function  of  Relative  Delay  and  Phase  of 

Multipath  [35]. 
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The  theoretical  basis  for  Figure  5  is  well  developed  in  [4],  The  nature  of  the  spread- 
spectrum  GPS  signal  causes  multipath  delayed  by  more  than  1.5  chips  to  be  rejected. 
However,  multipath  with  shorter  delays  may  cause  ranging  errors.  The  value  of  the  error 
due  to  multipath  can  be  either  positive  or  negative,  depending  on  the  phase  of  the 
multipath  relative  to  that  of  the  direct  signal.  The  amplitude  and  shape  of  the  error 
envelope  depends  on  the  relative  amplitude  of  the  multipath  signal.  As  the  amplitude  of 
the  multipath  increases,  the  error  bounds  in  Figure  5  grow  larger,  and  the  error  envelope 
becomes  increasingly  asymmetrical.  Therefore,  the  magnitude  of  the  ranging  error  due  to 
multipath  is  a  function  of  S,  a,  and^ .  The  spectral  characteristics  of  the  multipath 

error  are  dependent  on  <p  [28]. 

For  this  thesis,  the  focus  of  multipath  is  its  effect  on  the  airborne  platforms  and  the 
munition.  On  airborne  platforms,  the  relative  delay  ^is  limited  by  the  size  of  the 
airplane.  By  limiting  S ,  the  magnitude  of  the  multipath  ranging  error  is  also  limited.  A 
study  of  multipath  delay  on  Boeing  777,  747-400  and  737  aircraft  showed  a  20  to  60  cm 
1-a  ranging  error  using  pseudorange  measurements  [28]. 

To  determine  the  effects  of  multipath  on  a  munition,  a  study  on  multipath  effects  on  a 
satellite  is  used  as  a  point  of  comparison.  This  is  due  to  the  similarity  in  the  amount  of 
surface  area  in  which  a  signal  could  reflect  off  and  enters  the  antenna  as  delay.  A  study 
by  Godet  found  that  the  error  range  in  differential  phase  measurements  to  be  in  the 
millimeter  range  [15].  As  forementioned,  this  thesis  does  not  use  phase  measurements. 
However,  due  to  the  very  small  limit  on  the  amount  of  delay  possible  from  a  reflected 
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signal  on  a  munition,  code  multipath  is  expected  to  be  on  the  order  of  10  cm  or  less  for 
falling  munitions  [5]. 

Another  factor  in  munition  multipath  is  the  ground.  As  an  airplane  (or  munition) 
approaches  the  ground,  the  geometry  between  the  airplane,  the  satellite,  and  any 
multipath  sources  on  the  ground  is  changing  rapidly.  Hence  the  relative  phase,  cp ,  and 
the  relative  time  delay,  6 ,  are  likewise  changing  rapidly.  As  the  phase  continually 
changes,  the  multipath-induced  error  changes  rapidly  from  negative  to  positive,  making  it 
appear  noise-like  in  appearance  when  periodically  sampled  [28].  Smoothing  the  code 
measurements  with  the  carrier-phase  measurement  can  reduce  this  multipath. 

Multipath  is  generally  worse  for  fixed  receivers,  such  as  DGPS  ground  stations. 
Depending  on  the  ground  environment  (location  of  obstacles  from  which  signals  can 
reflect),  antenna  location,  and  receiver  type,  multipath  errors  could  be  anywhere  from  0 
to  5  meters  [5].  Since  this  research  does  not  use  a  ground  reference  station,  this  level  of 
multipath  error  is  not  a  concern. 

2.6.3.  Measurement  Noise  ( S  ) 

Measurement  or  receiver  noise  errors  are  due  to  the  measurement  processes  used 
within  the  receiver.  These  include  the  design  of  the  antenna,  the  method  used  for  the 
analog  to  digital  conversion,  the  correlation  processes,  and  the  tracking  loops  and 
bandwidths  [36].  Similar  to  the  multipath,  the  level  is  dependent  on  the  signal  that  is 
tracked.  The  C/A  code  measurement  noise  will  be  approximately  one  order  of  magnitude 
greater  than  P-code  measurements.  One-sigma  values  are  usually  on  the  order  of  1.5m 
for  the  C/A  code  and  20  cm  on  the  P-code  [19].  The  measurement  noise  due  to  carrier- 
phase  is  usually  between  1.2mm  and  1.6  mm  [30,31].  These  errors  are  uncorrelated,  and 
they  are  typically  modeled  as  pure  white  noise. 
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2.7.  DGPS  Effect  on  GPS  Errors 

Using  DGPS  is  known  to  produce  more  accurate  solutions  than  stand-alone  GPS. 
While  DGPS  may  reduce  or  remove  some  common  errors,  it  can  also  cause  some  errors 
to  be  added  or  even  amplified. 

Figure  6  shows  a  typical  DGPS  setup,  where  Measi  is  the  measurement  from  the 
satellite  to  the  reference  receiver  and  Meas2  is  the  measurement  from  the  satellite  to 
receiver  2. 


(Reference 

Receiver) 


Satellite 


Figure  6.  DGPS  Concept 


The  comparison  between  the  measurements  and  the  DGPS  effect  on  the  error  are  found  in 
Table  1[30]. 
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Table  1.  DGPS  Effect  on  GPS  Errors 


Error  Type 

Comparison  between  Measi  and  Meas2 

DGPS  Effect  on  Error 

Satellite  clock  error 

Same 

Removed 

Receiver  clock  error 

Different  (uncorrelated) 

Added 

Ephemeris  error^ 

Very  similar"' 

Reduced"' 

Ionospheric  delay 

Very  similar^ 

Reduced"' 

Tropospheric  delay 

Very  similar"' 

Reduced^ 

Multipath 

Different  (uncorrelated) 

Added  (and  amplified) 

Measurement  noise 

Different  (uncorrelated) 

Added  (and  amplified) 

^Effect  of  ephemeris  error  on  positioning  (actually  only  affects  the  calculated  range,  not 
the  actual  measurement) 

^Errors  grow  as  the  separation  distance  between  receivers  1  and  2  increases.  (The  errors 
are  the  same  and  are  removed  for  very  short  baselines  distances) 


The  satellite  clock  error  is  cancelled.  Differential  ephemeris,  ionospheric,  and 
tropospheric  errors  are  based  upon  the  distance  between  the  two  receivers.  At  a  close 
range,  the  differential  errors  are  reduced  because  of  the  correlation  between  the  errors, 
but  as  the  distance  grows,  the  amount  of  error  slowly  increases.  Multipath  and  the 
measurement  noise  are  uncorrelated  errors,  so  they  do  not  cancel.  In  fact  not  only  are  the 
errors  summed,  but  also  they  are  amplified  by  the  double  difference. 


2.8.  Synthetic  Aperture  Radar/Ground  Moving  Target  Indicator 

Synthetic  Aperture  Radar  (SAR)  produces  high-resolution  images  of  surface  areas 
and  has  the  ability  to  operate  in  all-weather  conditions.  SAR  can  be  used  on  many 
different  platforms,  including  airplanes,  the  Space  Shuttle,  and  satellites  [38].  Because 
radar  imagery  resolution  is  a  function  of  the  radar  sensor’s  aperture,  a  larger  aperture 
produces  higher  resolution  imagery.  A  SAR  uses  the  motion  of  the  airborne  platform  to 
synthesize  a  large  aperture  antenna  from  the  true,  smaller  aperture  antenna.  Typical  SAR 
sensors  provide  two  modes  of  operation:  search  and  spotlight.  In  search  mode,  a  SAR 
will  radiate  a  swath  of  land  providing  a  large  area  (and  usually  lower  resolution)  image. 
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In  spotlight  mode,  the  SAR  radiates  a  smaller  area  multiple  times  producing  a  higher 
resolution  image.  Figure  7  shows  a  SAR  mission  with  both  search  and  spotlight  mode. 
A  SAR  utilizes  typical  radar  techniques  by  measuring  the  time  between  the  transmission 
and  reception  of  a  SAR  signal  [23]. 


Spotlight  Mode 


Search  Mode 

Figure  7.  Synthetic  Aperture  Radar  Techniques[46] 


SAR  targets  are  typically  stationary.  However,  the  Ground  Moving  Target  Indicator 
(GMTI)  is  a  radar  mode  that  can  pickup  moving  targets.  The  radar  returns  echoes  from 
moving  vehicles  on  the  ground  are  separated  from  the  ground  clutter  on  the  basis  of  their 
Doppler  frequencies  [6]. 

2.8.1.  SAR  Errors 

For  this  research,  the  measurements  received  from  the  SAR  are  assumed  to  be  a  range 
and  range  rate.  In  early  AMSTE  studies,  SAR  measurements  were  range  and  range  rates 
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[6].  The  errors  for  the  range  can  occur  as  a  bias  and/or  a  noise,  where  biases  have  a  long 
correlation  time  and  noises  have  a  short  correlation  time. 

Some  bias-like  range  errors  include  range  clock,  and  atmospheric  refraction  [23]. 
Range  measurements  are  made  by  measuring  the  time  difference  between  the  transmitted 
pulse  and  its  return.  Consequently,  any  error  found  in  the  clock  rate  will  result  in  a  range 
measurement  error  [23].  In  general,  a  clock  runs  either  slightly  fast  or  slightly  slow  with 
very  slowly  changing  variations  in  the  rate  [23].  Given  the  assumption  of  no  time  latency 
in  this  thesis,  this  effect  goes  away.  The  atmospheric  refraction  error  is  similar  to  the 
troposphere  error  in  GPS.  It  is  a  function  of  how  the  particles  in  the  air  affect  the  path  of 
the  signal. 

Noises  on  the  range  can  include  quantization  errors  and  range  timing  errors.  The 
quantization  errors  result  from  resolution  errors.  The  range  timing  error  results  from  the 
time  delay  needed  to  perform  SAR  processing.  Both  of  these  errors  are  modeled  as  white 
noise. 

The  range  rate  errors  are  primarily  noise-like  in  character.  The  radar’s  receiver  noise 
and  internal  motion  resolution  errors  can  add  various  levels  of  error.  Calibration,  if  it  is 
not  precise,  can  add  a  bias-like  error  to  the  range  rate.  More  quantitative  details  of  SAR 
errors  are  given  in  Chapter  3. 

2.9.  System  Kalman  Filter 

The  Kalman  filter  is  an  optimal  recursive  data  processing  algorithm  [25,26].  The 
filter  is  optimal  referring  that  the  entire  history  of  measurements  available  to  the  filter  is 
processed  and  incorporated  as  information  in  the  form  of  dynamics  model,  measurement 
model,  and  statistical  descriptions  of  noises  and  uncertainties.  The  filter  is  recursive  in 
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that  the  filter  does  not  require  all  of  the  previous  data  to  be  kept  in  storage  and 
reprocessed  every  time  a  new  measurement  becomes  available.  The  Kalman  filter 
processes  all  available  measurements  of  interest,  regardless  of  their  accuracy,  based  on 
knowledge  of  the  system  and  measurement  dynamics,  the  statistical  description  of  the 
system  noises,  measurement  errors,  and  model  uncertainties  [25,26].  The  models  may  be 
linear,  in  certain  cases,  but  most  system  models  are  nonlinear  in  nature.  For  nonlinear 
system  models,  the  extended  Kalman  filter  (EKF)  is  implemented.  The  extended  Kalman 
Filter  is  not  “optimal,”  but  only  a  first-order  approximation  to  the  (infeasible,  infinite¬ 
dimensional)  optimal  non-linear  filter  [25,26].  The  EKF  linearizes  the  nonlinear  system 
model  about  a  nominal  point  in  the  trajectory,  utilizing  the  linear  assumptions  and 
equations  described  in  the  following  sections.  During  operation,  an  EKF  is  relinearized 
based  on  the  most  current  optimal  estimate  of  the  variable  of  interest. 

2.9.1.  State  and  Measurement  Model  Equations 

In  developing  the  Extended  Kalman  filter  for  this  research,  let  the  system  model  be 
defined  as  a  state  equation  in  the  following  form  in  Equation  (2-12), 

x(f,.)  =  a>(t,.,Vi)x(t,._i)-HGa(i,._,)Wd(0  (2-12) 

where 

x(r . )  =  state  vector  at  time  q 

x(q._, )  =  state  vector  at  time  tj.i 

^(q  .^-1 )  =  linear  discrete  state  dynamics  matrix 

G  j  (q_i )  =  noise  distribution  matrix  (the  identity  matrix  for  this  system) 

w  j  (f)  =  white  Gaussian  noise  with  a  mean  value  of  zero  and  strength 


E[w^(t)]  =  0 


(2-13) 


E[w^(t-)w/ (tj)]  = 


{QAti)for  t,=tj 

jo  for 


(2-14) 
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where 

Qd  (h )  =  covariance  of  w ^  (?,. ) 

The  discrete  time  measurements,  z(t,  ),  are  modeled  in  Equation  (2-15) 

z(t,)  =  h[x(?,),t,]-hv(f,)  (2-15) 

where 

h[x(?; ),  ?,  ]  =  nonlinear  function  of  the  state  vector  and  time 
v(?, )  =  discrete  time  measurement  noise  vector  with  a  zero-mean  white  noise 

process,  which  is  independent  of  w(t,.) ,  and  having  covariance  R(t,  ) 
defined  by: 

£[v(0]  =  0  (2-16) 


^[v(r,)v^(t,)] 


for  ti=tj 
jo  for  r,.  lit  tj 


(2-17) 


2.9.2.  Measurement  Model  Linearization 

Since  the  measurements  are  nonlinear,  the  EKF  filter  equations  must  be  linearized  to 
produce  the  prediction  of  the  measurement  vector  z(?,.)  before  it  arrives,  to  process  a 
measurement  update  cycle  of  the  filter.  The  linearization  is  performed  by  taking  the 
partial  derivative  of  h[x(r, ),/,.]  to  produce  H[x„  (/,.),?,.] 


H[t.;x(tr)]  = 


^[x(0,t]| 


(2-18) 


lx=x(;f ) 


Normally,  if  the  true  state  trajectory  differs  from  the  desired  state  trajectory,  large 
errors  could  occur.  The  EKF  reduces  this  effect  by  allowing  a  relinearization  about  the 
most  recent  state  estimate  as  shown  in  Equation  (2-18),  as  opposed  to  the  linearized 
Kalman  filter  which  only  uses  the  nominal  state  value.  Using  the  EKF  method  allows  for 
the  declaration  of  a  new  nominal  pre-computed  trajectory  (emanating  from  the  most 
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recent  state  estimate)  at  every  estimate.  This  ensures  deviations  from  the  nominal 
trajectory  remains  small  as  long  as  the  error  model  stays  accurate. 

2.9.3.  Extended  Kalman  Filter  Propagate  and  Update  Equations 

This  section  addresses  the  EKF  propagate  and  update  equations  implemented  in  this 
research  effort.  Equation  (2-12)  is  the  discrete-time  dynamics  model  and  Equation  (2-15) 
is  the  discrete-time  measurement  model.  As  previously  mentioned,  a  Kalman  filter  is  a 
recursive  algorithm.  There  are  two  steps  involved  in  this  recursion:  propagation  and 
update.  The  state  estimate,  x(0 ,  and  the  covariance  of  that  estimate,  P(0 ,  are  both 
propagated  from  the  last  time  sample,  f,.;,  and  updated  at  every  time,  .  Sampled  data 
EKF  equations  utilize  the  following  notations: 

tr  -  value  of  a  variable  after  propagation  from  but  just  prior  to  a  measurement 

update  at  time  t,-. 

-  value  of  a  variable  after  propagation  from  ti.i  and  the  measurement  update  at  time 
t,. 

The  subscript  i  is  used  to  describe  the  discrete  time  points  when  measurements  are 
available.  Using  these  time  notations,  the  state  estimates  x(t”)  and  covariance  values 
P(?,)  are  propagated  from  to  using  the  following  discrete  equations: 

x(r-)  =  a)(t,-t,.,)x(C,)  (2-19) 

Pit; )  =  a)(t,  - )P(C,  )o"  (t,  -  ) + Q,  (2-20) 

These  equations  are  for  time-invariant  system  model  and  stationary  noise  model.  When 
discrete  time  measurements,  z,.,  become  available,  the  EKF  update  cycle  is  performed 
using  the  following  equations: 
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K(f,0  =  P(rr)H^[?,;x(rr)]{H[?,;x(f-)]Pa-)H^[f,;x(^^^ 


(2-21) 


where 

K{t,)  =  Kalman  Filter  Gain 


x(t;)  =  x(t:)  +  K(i,){z,  -h[x(f-),f,]}  (2-22) 

P(t;)  =  (I-K(r,OH[t,;x(t:)])P(i-)(I-K(i,)H[i,.;x(t:)])"  +K(t,.)R(t,)K(r,)"  (2-23) 

For  more  details  and  a  complete  derivation  of  the  above  equations,  the  reader  is  referred 
to  [25,26]. 

2.10.  Environment  Model 

The  Environment  Model  will  be  the  truth  model  for  the  system  model.  A  truth 
model  results  from  a  detailed  analysis  of  the  system,  and  it  includes  errors  and  other 
characteristics  that  may  be  assumed  negligible  by  the  system  model.  It  should  also  be  the 
best  possible  representation  of  the  “real  world.”  Typically,  a  truth  model  is  contains 
many  more  state  variables  than  the  system  filter  design  model  and  a  higher  level  of 
accuracy.  However  due  to  the  nonlinear  deterministic  nature  of  some  errors  (troposphere 
in  particular),  meaning  that  they  can  be  found  with  known  parameters  and  are  not  truly 
random  in  nature,  the  “truth”  environment  model  is  not  only  driven  by  white  noise.  In 
addition,  the  environment  model  will  begin  with  profiles  of  exact  position,  velocity  and 
acceleration  (if  necessary)  for  each  platform,  the  munition,  satellites  and  target.  Random 
values  representing  the  errors  in  the  measurements  will  be  added.  For  example,  the  true 
range  between  satellite  #1  and  sensor  platform  #1  at  a  given  time  is  known  from  the 
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profiles.  GPS  errors  such  as  satellite  clock  and  ephemeris  along  with  atmospheric 
(troposphere,  ionosphere)  and  receiver  errors  will  be  added,  producing  a  known 
pseudorange.  Other  neglected  errors  (neglected  in  the  filter  design  model)  such  as 
multipath  and  the  unmodeled  troposphere  errors  will  also  be  added.  This  will  show  what 
effect  not  modeling  them  in  the  filter  design  model  has  on  the  resulting  total  error.  The 
actual  elements  of  the  environment  model  will  be  discussed  in  detail  in  the  next  chapter. 
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3.  Modeling  Methodology 


3.1.  Overview 

This  chapter  focuses  on  the  development  of  the  true  environment  model  and  the 
system  models,  integration  methods,  and  simulation  techniques  used  in  this  research 
effort.  Section  3.2  describes  the  true  environment  and  all  the  models  that  produce  the 
GPS  measurements  and  the  SAR  range  and  range  rate  measurements  used  by  the  Kalman 
filter.  Section  3.3  examines  the  elements  of  the  Kalman  filter.  Section  3.4  provides  the 
equations  and  calculations  necessary  to  examine  to  error  in  the  range  from  the  bomb  to 
the  target. 

3.2.  True  Environment  Modeling 

The  true  environment  model  is  illustrated  in  Figure  8: 


GPS 


Figure  8.  True  Environment  Model 
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The  true  environment  consists  of  two  sensor  platforms  (a  Global  Hawk  and  a 
JSTARS),  the  munition  delivered  by  a  fighter,  a  target,  and  the  GPS  constellation  of 
satellites  in  view  to  the  receivers  on  the  sensor  platforms  and  munition.  The  lightning-like 
lines  represent  the  pseudorange  measurements  from  the  GPS  satellites  to  the  GPS 
receivers.  The  gray  dashed  lines  from  the  sensors  to  the  bomb  and  target  represent  the 
SAR  range  and  range  rates.  The  dotted  line  from  the  fighter  to  the  bomb  represents  the 
flight  path  of  the  bomb  to  impact  the  target.  The  three-way  axis  represents  the  local- 
level,  East-North-Up  (ENU)  coordinate  frame.  The  line  from  the  bomb  to  the  target 
represents  the  range  between  the  two  and  is  the  primary  focus  of  the  simulations:  to 
minimize  the  error  in  the  range  from  the  bomb  to  the  target.  The  scenario  is  750  seconds 
in  length  and  was  based  on  preliminary  work  for  the  AMSTE  program  [17,43]. 

3.2.1.  Sensors  Environment  Models 

The  models  of  the  two  sensors,  also  referred  to  as  Sensor  1  and  2,  were  built  in 
MATLAB  using  PROGEN,  a  flight  profile  generator  [2].  PROGEN  is  a  local-level  flight 
trajectory  generator  suitable  for  short-distance,  short-duration  flights.  Sensor  1  is  at  an 
altitude  of  10  kilometers,  and  sensor  2  at  16.7  kilometers,  and  both  fly  the  same 
rectangular  racetrack  shape  trajectory.  Table  2  defines  this  shape.  They  are  a  separated 
by  a  nearly  constant  look  angle  to  the  target  of  40  degrees.  Figure  9  shows  a  top  view  of 
the  two  sensor  flight  profiles  along  with  the  top  view  of  the  bomb  and  target.  The  target 
is  offset  from  the  profile  to  allow  it  to  be  seen. 
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Table  2.  Racetrack  shape  of  Sensor  1  and  2 


Fly  straight 

180  sec 

Roll  right 

45  deg 

Turn  90 

deg 

Roll  back  to  level 

Travel  for 

90  sec 

Roll  right 

45  deg 

Turn  90 

deg 

Roll  back  to  level 

Travel  for 

180  sec 

Roll  right 

45  deg 

Turn  90 

deg 

Roll  back  to  level 

Travel  for 

90  sec 

East  North  view  of  Sensors, Bomb  and  Target 


Figure  9.  East-North  view  of  Sensor  1  and  Sensor  2  Flight  Profiles 


3.2.2.  Munition  Environment  Model 

The  profile  for  the  munition,  or  bomb,  is  described  as  a  bomb  being  carried  by  a 
fighter  that  flies  at  a  constant  altitude  of  7620  meters  for  630  seconds.  The  bomb  is  then 
dropped  such  that  it  impacts  the  target  from  a  nearly  vertical  direction.  By  coming  in 
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vertically,  the  error  in  the  Up  direction  of  the  target,  which  will  be  discussed  later, 
becomes  less  of  a  factor  in  accurately  destroying  the  target.  Figure  10  gives  the  profile  of 
the  bomb  altitude  for  the  scenario. 


Bomb  Altitude  Projection 


Figure  10.  True  Environment  Bomb  Altitude 

3.2.3.  Target  Environment  Model 

The  target  model  is  very  simple.  To  maintain  the  assumptions  of  only  one  visible  and 
constantly  moving  vehicle,  the  target  begins  at  the  origin  of  the  ENU  coordinates  and 
moves  eastward  at  a  constant  velocity  of  24.2  m/s.  The  trajectories  are  designed  so  that 
the  bomb  impacts  the  target  at  t=750  seconds.  This  is  important  because  this  research  is 
only  looking  at  the  navigation,  not  the  guidance  of  the  bomb. 
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3.2.4.  GPS  Environment  Model 


The  GPS  environment  model  is  built  by  a  randomly  picked  ephemeris  data  file.  The 
ephemeris  data  file  is  randomly  chosen  to  provide  variety  in  the  satellite  constellations 
for  Monte-Carlo  simulations.  The  ephemeris  file  is  used  to  calculate  the  positions  of  the 
satellites  that  are  visible  to  the  three  receivers.  The  next  few  sections  explain  the  GPS 
and  SAR  environment  true  ranges  and  range  rates  and  the  errors  that  are  added  to  the  true 
ranges/range  rates  to  produce  the  measurements. 

3.2.4.I.  Environment  Model  True  Ranges  and  Range  Rates 

This  section  shows  the  equations  used  to  build  the  true  GPS  ranges  and  the  SAR 
range  and  range  rates.  The  GPS  ranges  are  found  by  calculating  the  distance  between 
the  GPS  satellites  and  the  receiver  positions  [30].  This  is  repeated  for  each  receiver  and 
the  satellites  visible  to  the  receiver. 

GPS pseudoran,e  =  +  ()'.v  "  )'  +  +  cS  t,  -  cS  t,, 

+  +  cS  +  cS  t„p  -h  S  r^p^ 


where 


^SV  ^  y  SV  ^  ^5V 

=  EMU  positions  of  the  satellite 

=  ENU  positions  of  the  receiver 

cSt^ 

=  receiver  clock  error 

cSt^^ 

=  satellite  clock  error 

Cp, 

=  error  due  to  broadcast  ephemeris 

^  trop 

=  delay  due  to  the  troposphere 

O 

a 

o 

=  delay  due  to  the  ionosphere 

noise 

=  receiver  noise 

=  multipath  error 

=  ephemeris  error 

The  model  for  each  error  will  be  discussed  in  the  next  sections. 
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3.2.4.2.  GPS  Troposphere  Environment  Model 

As  mentioned  in  chapter  2,  the  modified  Hopfield  model  was  used  to  determine  the 
tropospheric  error.  The  modified  Hopfield  model  uses  the  receiver’s  height,  elevation, 
and  relative  humidity  to  estimate  the  error  that  needs  to  be  subtracted  from  the 
pseudorange.  There  is  usually  a  residual  error  after  the  model  has  been  applied.  This 
residual  is  the  error  that  gets  applied  to  the  true  range  in  the  research.  To  build  the 
residual  error,  a  true  relative  humidity  between  0%  and  100%  is  randomly  selected  and 
the  true  tropospheric  error  is  determined.  The  model  is  run  again,  but  this  time  the 
relative  humidity  is  set  at  50%.  This  value  is  chosen  because  normally  a  value  of  50%  is 
used  when  the  true  humidity  is  not  known.  The  two  model  estimates  are  then  subtracted 
from  each  other  to  represent  the  residual  error.  This  error  can  be  10"^  to  4x10  ^  meters  of 
error  for  the  sensors  and  0.1  to  0.6  meters  of  error  for  the  bomb.  The  amount  of  error  is 
dependent  on  the  height  of  the  receiver. 

3.2.43.  GPS  Environment  Satellite  Clock  Model 

The  correction  coefficients  for  the  satellite  clock  error  are  transmitted  as  part  of  the 
navigation  message  (afo,afi,af2).  Once  the  corrections  have  been  applied  in  the  real  world, 
the  resulting  1-a  residual  error  is  typically  10ns  or  3m  [19].  In  the  environment  model, 
the  satellite  clock  error  is  modeled  as  a  constant  bias  for  the  duration  of  the  scenario.  For 
the  length  of  time  of  the  scenario  (750  seconds),  this  error  does  not  change  significantly. 
The  bias,  chosen  independently  for  each  satellite,  is  a  normally  distributed  random 
variable  with  a  standard  deviation  of  10ns.  This  error  is  multiplied  by  the  speed  of  light 
to  convert  to  meters. 
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3.2.4.4.  GPS  Environment  Ephemeris  Model 

The  satellite  ephemeris  error  model  uses  the  satellite  positions  sent  in  the  navigation 
message  and  precise  satellite  positions  from  post-processed  data.  The  difference  between 
the  transmitted  navigation  message  position  and  precise  positions  is  the  ephemeris  error. 
These  errors  are  transformed  from  the  downtrack,  crosstrack  and  radial  directions  to  the 
ENU  coordinate  frame.  Calculating  the  projection  of  the  error  onto  the  unit  line  of  sight 
of  the  receiver  gives  the  appropriate  level  of  error,  about  3m  (la)  [38]. 

3.2.4.5.  GPS  Environment  Receiver  Clock  Model 

The  environment  receiver  clock  model  is  described  in  Figure  11.  The  strength  of  the 
white  noises,  uv  and  us,  are  functions  of  Allan  variance  parameters  [31].  Their  values  are 
calculated  in  the  following  equations: 

Sf=2xh^  (3-2) 

S^=S7r^xh_^  (3-3) 

where 

Sj-  =  Strength  of  white  noise  u7 

Sg  =  Strength  of  white  noise  u8 

K  ’  ^-2  =  Allan  variance  parameters  typical  quartz  crystal,  Hq  =  2E-19,  h_2  =  2E-20 

White  noise  u 


Figure  11.  Environment  Receiver  Clock  model 
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3.2.4.6.  GPS  Environment  Receiver  Noise  Model 


The  noise  of  a  receiver  is  uncorrelated  with  the  noise  of  any  other  receiver  [31].  Each 
receiver’s  noise  is  modeled  as  a  white  Gaussian  noise  with  a  standard  deviation  of  0.1m, 
which  is  a  noise  value  typical  for  P-code  measurements  [19]. 

3.2.4. 7.  GPS  Environment  Multipath  Model 

The  multipath  environment  model  for  both  the  sensor  platforms  and  the  bomb  is 
modeled  as  a  first  order  Gauss-Markov  process.  This  model  is  consistent  with  [5].  The 
first-order  Markov  model  parameters  are  a  standard  deviation  1  meter  with  a  time 
constant  of  30  seconds.  The  bomb  has  multipath  of  10cm  (la)  and  a  time  constant  of  30 
seconds.  Figure  12  gives  an  example  of  multipath  for  an  airplane. 


Multipath  for  Sensor  1 


Figure  12.  Multipath  Error  Example  for  Sensor  1 
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3.2.5.  SAR/GMTI  Environment  Range  and  Range  Rate  Models 

The  SAR/GMTI  trae  ranges  to  the  bomb  and  target  are  calculated  from  the  positions 
of  the  profiles.  The  SAR  true  range  is  calculated  in  the  same  manner  as  the  GPS  range, 
except  it  is  the  difference  between  the  sensor  and  the  bomb  or  the  target  [23,46].  The 
errors  added  to  the  range  and  range  rate  to  produce  the  SAR  measurements  are  discussed 
in  the  following  sections. 

=  ^(xp-xy+iyp-yf+izp-z,f  +  5  trop  +  (3-4) 

where 


>  y  p  ’  ^p 

=  ENU  positions  of  the  sensor 

x,,y,,z, 

=  ENU  positions  of  the  target/bomb 

S  trop 

=  SAR  atmospheric  error 

V 

noise 

=  SAR  range  noise 

The  SAR  range  rate  is  a  function  of  the  difference  of  velocities  projected  along  the  line  of 
sight 


SAR 


P  -P 

Ip  X, 


rangerate 


^(Xp  -x^  +{yp-yy  +  (zp  -z,y 


'(Vp-v,  )  + 


(3-5) 


where 

Pp  =  position  vector  of  the  sensor 
Pj  =  position  vector  of  the  target/bomb 
\p  =  velocity  vector  of  the  sensor 
Vj  =  velocity  vector  of  the  target/bomb 
Vnoise=  SAR  range  rate  noise 
•  =  dot  product 


3.2.5.I.  SAR/GMTI  Environment  Atmosphere  Range  Error 

The  environment  range  atmosphere  error  is  similar  to  the  tropospheric  error  for  the 
GPS  signal.  In  the  case  of  the  SAR  range,  the  troposphere  not  only  delays  the 
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transmission  of  the  signal,  but  also  the  return.  The  error  is  modeled  in  the  same  manner 
as  the  troposphere  error  is  modeled,  however  an  additional  0.1”  error  (la)  is  added  to  the 
elevation  angle  to  model  the  mapping  function  errors.  The  mapping  function  errors  are 
due  to  the  interpolation  of  the  model  for  elevation  angle.  The  difference  of  the  two 
calculated  atmosphere  errors,  which  is  multiplied  by  2  to  account  for  the  transmission  and 
return,  becomes  the  error  due  to  the  atmosphere.  This  is  illustrated  in  Figure  13. 

The  process  is  repeated  for  the  bomb.  The  bomb  starts  at  a  higher  altitude,  where 
there  are  lower  atmospheric  errors.  As  the  bomb  drops  altitude  and  approaches  the  target, 
the  value  of  the  bomb  and  target  atmospheric  errors  converges.  Figure  14  gives  a 
graphical  depiction  of  a  sample  case.  The  changes  in  the  magnitude  of  error  are  due  to  the 
changes  in  elevation  angle. 


Figure  13.  SAR  Atmosphere  error  model  (value  is  doubled  for  transmission  & 

reception) 
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SAR/GMTI  Environment  Atmosphere  Error 


Figure  14.  SAR/GMTI  Environment  Atmosphere  Error 


3.2.5.2.  SAR/GMTI  Environment  Range/Range  Rate  Noise  Error 

There  are  a  number  of  factors  that  add  noise-like  errors  to  a  range  or  range  rate, 
including  timing  and  quantization  errors  [23].  The  modeling  of  the  range  and  range  rate 
noises  come  from  the  work  done  in  [23,  46].  Taking  the  root  mean  square  of  all  the  noise 
terms,  the  range  noise  error  is  3.17  meters  (Ic)  and  the  range  rate  noise  error  is  0.04 
meters/sec  (la). 

3.2.6.  True  Environment  Summary 

All  of  the  components  and  error  models  have  been  described.  Table  3  summarizes 
the  true  ranges  and  the  associated  errors  added  to  the  ranges  to  produce  the  measurements 
used  by  the  System  Model. 
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Table  4.  Filter  Model  States 


State 

# 

Symbol 

Description 

State 

# 

Symbol 

Description 

1 

East  position  of  Sensor  1  (m) 

22 

East  velocity  of  Bomb  (m/s) 

2 

ypi 

North  position  of  Sensor  1  (m) 

23 

ys 

North  velocity  of  Bomb  (m/s) 

3 

Up  position  of  Sensor  1  (m) 

24 

Zb 

Up  velocity  of  Bomb  (m/s) 

II 

Xp^ 

East  velocity  of  Sensor  1  (m/s) 

25 

Xq 

East  acceleration  of  Bomb 
(m/s^) 

H 

North  velocity  of  Sensor  1 
(m/s) 

26 

yB 

North  acceleration  of  Bomb 
(m/s^) 

6 

Zpi 

Up  velocity  of  Sensor  1  (m/s) 

27 

Zb 

Up  acceleration  of  Bomb  (m/s"^) 

7 

Xpi 

East  acceleration  of  Sensor  1 
(m/s^) 

Xj 

East  position  of  Target  (m) 

8 

North  acceleration  of  Sensor  1 
(m/s^) 

yp 

North  position  of  Target  (m) 

9 

Zpi 

Up  acceleration  of  Sensor  1 
(r^s^) 

Xj 

East  velocity  of  Target  (m/s) 

10 

^P2 

East  position  of  Sensor  2  (m) 

31 

% 

North  velocity  of  Target  (m/s) 

11 

>'P2 

32 

Xj 

East  acceleration  of  Target 
(m/s^) 

12 

^P2 

Up  position  of  Sensor  2  (m) 

33 

yr 

North  acceleration  of  Target 
(m/s^) 

wm 

Xp2 

East  velocity  of  Sensor  2  (m/s) 

34 

Clock  bias  Sensor  1  (m) 

m 

North  velocity  of  Sensor  2 
(m/s) 

35 

Clock  drift  Sensor  1  (m/s) 

15 

Up  velocity  of  Sensor  2  (m/s) 

36 

Clock  bias  Sensor  2  (m) 

H 

Xp2 

East  acceleration  of  Sensor  2 
(m/s^) 

37 

cSi^ 

Clock  drift  Sensor  2  (m/s) 

■ 

ypi 

North  acceleration  of  Sensor  2 
(m/s^) 

38 

cJtpj 

Clock  bias  Bomb  (m) 

H 

‘^P2 

Up  acceleration  of  Sensor  2 
(r^s^) 

39 

C<5tp2 

Clock  drift  Bomb  (m/s) 

19 

Xq 

East  position  of  Bomb  (m) 

40 

sfpi 

SAR  scale  factor  of  target  - 
platform  1  (m/m) 

20 

yp 

North  position  of  Bomb  (m) 

41 

P2 

SAR  scale  factor  of  target  - 
platform  2  (m/m) 

21 

Zb 

Up  position  of  Bomb  (m) 

42+ 

Satellite  #_  errors  (m) 

The  next  sections  describe  the  dynamics  models  for  these  states  and  give  the  initial 
value  and  covariances.  The  dynamics  models  for  these  states  can  be  in  either  continuous 
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form  or  an  equivalent  discrete-time  form.  Both  forms  are  shown  in  the  following 
equations. 


X  =  Fx  +  Bu  +  Gw 

x(r,. )  =  Ox(t,.., )  +  B^uCt;.! )  +  Wj  ) 


(3-63-7) 


where 

F,  <I>  =  continuous,  discrete-time  forms  of  the  dynamics  matrix 

B,Bd  =  continuous,  discrete-time  forms  of  the  input  matrix  (zero  matrix  for  study) 
=  continuous,  discrete-time  forms  of  the  vector  of  input  states(zero  for  this 
study) 

w,Wd(t|.;)  =  continuous,  discrete-time  forms  of  the  vector  of  driving  white  noise 
For  these  descriptions,  the  state  vector  will  be  split  into  segments.  For  example,  the 
position,  velocity,  and  acceleration  states  of  sensor  1  are  one  segment.  The  ^  and  the  Qd 
matrices  are  shown  below  in  block  form. 


o= 


FOGMAp, 

0 

0 

0 

0 

0 

0 


0 

FOGMAp^ 

0 

0 

0 

0 

0 


0 

0 

FOGMAp 

0 

0 

0 

0 


0 

0 

0 

FOGMAp 

0 

0 

0 


0  0 

0  0 

0  0 

0  0 

Clock  0 

0  SAR^j 

0  0 


0 

0 

0 

0 

0 

0 

SAF 


Err. 


(3-8) 


where 

FOGMAp^ ,  FOGMAp^ ,  FOGMAp  =  9-by-9  3-D  First-Order  Gauss-Markov 

Acceleration  models  for  the  sensors,  bomb 
FOGMAp  =  6-by-6  2-D  First-Order  Gauss-Markov  Acceleration  model  for  the 
target 

Clock  =  three  2-by-2  Clock  bias/Clock  drift  models  for  the  receiver  clocks 

=  two  SAR  scale  factors  for  atmosphere  error  from  sensor  1/2  to  target 

SATp^^  =  N-by-N  matrix  of  satellite  error  models  (N=  #  of  satellites) 
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0 


Q 


dpi 


0 

0 

0 

0 

0 

0 


0 

QrfP2  ^ 

0  Q.. 

0  0 

0  0 

0  0 

0  0 


0  0 

0  0 

0  0 

Qrf7  0 

^  ^  d  clock 

0  0 

0  0 


0  0 

0  0 

0  0 

0  0 

0  0 

dSAR  ^ 

0  Q-s,. 


(3-9) 


where 

Q.P).Q.P2^Q.B=9-by  -9  Matrices  for  the  white  noise  on  the  acceleration  states 
=  6-by-6  Matrix  for  the  white  noise  on  the  acceleration  states 
Q^ciock  ~  6-by-6  Matrix  for  the  white  noise  of  the  receiver  clock  bias/clock  drift 
states 

QdsAR  -  2-by-2  Matrix  for  the  white  noise  of  the  SAR  modeled  random  walk 

=  N-by-N  matrix  for  the  white  noise  of  the  Gauss-Markov  satellite  models 


3.3.2.  Filter  Dynamics  Model  of  Sensor  1,  Sensor  2,  and  the  Bomb 

The  dynamics  model  for  each  of  the  two  sensors  and  the  bomb  is  a  9-by-9  3-D 
FOGMA.  The  position  states  are  an  integral  of  the  velocity  states,  the  velocity  states  are 
an  integral  of  the  acceleration  states  and  the  acceleration  states  are  each  modeled  as  a 
first-order  Gauss-Markov  process.  The  continuous-time  form  of  the  FOGMA  for  a 
sensor  or  bomb  is  represented  by  Fi. 
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0  0  0  1  0  0 
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(3-10) 


where 

T  =  time  constant  for  Gauss-Markov  process 

For  sensor  1  and  sensor  2  the  time  constant  is  3  seconds.  For  the  bomb  it  is  5 
seconds.  These  values  are  based  upon  the  anticipated  correlation  between  the 
acceleration  at  the  current  time  applied  to  the  next  time  [25,26,311- 

In  this  research,  the  model  is  implemented  in  discrete  time.  There  are  many  ways  to 
calculate  the  0  matrix  from  the  Fi  matrix.  Because  Fi  is  non-time  varying,  the  matrix 
exponential  method  is  used. 


^  t.  - 1._, 


(3-11) 

(3-12) 


The  characteristics  of  the  white  Gaussian  noise  values  for  the  FOGMA  are 
represented  first  in  continuous-time  form  with  the  following  equations: 

F[w]  =  0  (3-13) 
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E{w{t)v/  (?  +  t)}  =  it)  = 


00000000 

00000000 

00000000 

00000000 

00000000 

00000000 

000000^0 

0000000^ 

00000000 


d{T) 


(3-14) 


2(7^ 

~T~ 


(3-15) 


where 

q  =  strength  of  the  white  noise  on  the  acceleration  state 

2 

cr  =  standard  deviation  of  acceleration  in  meters/sec 
T  =  time  constant 

The  cr  value  for  sensors  1  and  2  is  8  meters  and  the  value  is  15  meters  for  the  bomb. 


As  with  the  dynamics  matrix,  the  covariance  matrix  needs  to  be  transformed  to 


discrete  time.  This  is  done  with  the  following  equation: 

Q,  =Jo(x)Qa»"(T)Jx 

0 


(3-16) 
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(3-17) 
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where 


The  initial  values  of  the  position  states  for  both  sensors  and  the  bomb  are  the  true 
positions  with  9m  (la)  of  error  added.  The  initial  values  of  the  velocity  and  acceleration 
are  zero.  Ppo  describes  the  initial  covariance 

'lOO^  0  0  0  0  oooo' 

0  100"  0  0  0  0  0  0  0 

0  0  100"  0  0  0  0  0  0 

0  0  0  400"  0  0  0  0  0 

P^O  =  0  0  0  0  400"  0  0  0  0  (3-18) 

0  0  0  0  0  400"  0  0  0 

0  0  0  0  0  0  20"00 

0  0  0  0  0  00  20"0 

00  0  0  0  000  20" 

The  large  initial  values  insure  that  the  Kalman  filter  will  place  more  weight  on  the  initial 
measurements  rather  than  the  unknown  initial  state  [31]. 

3. 3.3.  Filter  Dynamics  Model  for  Target 

The  dynamics  model  for  the  target  is  also  a  FOGMA  model  except  that  it  is  only  a  6- 
by-6  matrix  for  a  2-D  case  (East  and  North).  A  two-state  model  is  being  implemented 
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because  there  are  only  two  position  measurements  available  (sensor  1  and  sensor  2 
ranges).  The  same  matrices  and  equations  apply  as  above  minus  the  z  states  for  position, 
velocity  and  acceleration.  The  time  constant,  or  T ,  for  the  target  is  6  seconds  and  the  cr 
value  is  5  meters/sec^.  The  initial  state  values  and  covariances  are  the  same  as  that  for 
the  sensors  and  bomb. 

The  Up  direction  of  the  target  is  modeled  as  a  Digital  Terrain  Elevation  Data  (DTED) 
measurement.  DTED  is  an  estimate  of  the  height  of  a  position  given  the  longitude  and 
latitude  coordinates.  With  the  length  of  the  scenario  short,  the  error  in  height  is 
represented  by  constant  random  bias  with  a  1-g  error  of  10m.  For  the  target,  the  velocity 
and  acceleration  in  the  Up  direction  is  assumed  to  be  zero. 

3.3.4.  Filter  Model  for  Sensors  1  &  2  and  Bomb  Receiver  Clocks 

The  filter  model  for  the  GPS  receiver  clock  is  the  same  model  as  that  used  in  the  true 
environment  model  (as  shown  in  Figure  11).  The  system  model  is  a  two-state  model  of 
the  clock  bias  and  clock  drift.  The  clock  bias  accounts  for  a  large,  non-white  error  in  the 
code  measurements  [31].  While  the  clock  drift  is  not  needed  for  measurement 
incorporation,  is  it  useful  for  propagating  the  clock  bias  forward  in  time  and  preventing  a 
consistent  bias  in  all  the  code  measurements  [31].  Equations  (3-19)  and  (3-20)  represent 
the  continuous-time  form  of  the  dynamics  model  that  is  transformed  to  the  discrete  model 
in  the  same  manner  as  above.  This  model  is  repeated  for  each  receiver. 


“  clock 
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where 


Qci  =  0-036 


m 


qc2  =  0.141 


m 


The  initial  value  for  all  clock  bias  and  clock  drifts  are  zero  and  the  initial  covariance  is 
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0  1000" 


(3-21) 


3.3.5.  System  Model  for  SAR  scale  factor  for  the  target 

The  System  model  for  the  SAR  range  atmosphere  error  to  the  target  incorporates  a 
scale  factor  error.  This  model  was  implemented  in  [23,46].  Initially  the  scale  factor  was 
applied  to  both  the  target  and  bomb.  However,  after  some  initial  runs,  the  estimation  of 
the  atmosphere  error  was  significantly  off,  because  there  was  very  little  atmospheric  error 
between  the  sensor  and  the  bomb.  By  redefining  the  state  to  be  just  the  sensor-to-target 
atmospheric  error,  the  results  gave  more  reasonable  errors.  As  a  result,  the  measurement 
model  for  the  sensor-to-bomb  SAR  range  measurements  did  not  include  the  atmospheric 
error  term  until  the  bomb  was  in  the  vicinity  of  the  target  (defined  as  below  4000m 
altitude.)  This  is  not  the  best  model,  and  associated  concerns  will  be  covered  in  Chapter 
5-Conclusion  and  Recommendations. 

The  scale  factor  is  modeled  as  random  walk,  as  shown  in  equations  (3-22)  and  (3-23). 


F  = 

*  SAR 


0  0 
0  0 


(3-22) 


Qs/i/?  ~ 


1x10'®  0 

0  1x10'^ 


(3-23) 


6  2 

The  initial  value  is  zero  and  the  initial  covariance  is  (30x  10' )  [23,46]. 
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3.3.6.  Filter  Model  for  Common  GPS  Errors 

Here,  the  filter  models  errors  in  the  GPS  measurements  that  are  common  to  all  of  the 
receivers:  satellite  clock  and  satellite  ephemeris.  By  estimating  the  common  errors  for 
the  satellites  and  the  receiver  clock,  the  differential  calculation  is  performed  implicitly 
within  the  filter.  Not  estimating  the  errors  provides  a  stand-alone  GPS  solution.  In  the 
stand-alone  GPS  solution  case,  the  satellite  clock  and  ephemeris  errors  are  not  cancelled 
or  reduced,  but  left  as  part  of  the  pseudorange  measurement.  This  fact  will  be  shown  in 
the  next  chapter.  The  number  of  states  is  dependent  on  the  number  of  total  satellites 
visible  to  the  receivers.  The  common  error  is  modeled  as  a  Gauss-Markov  process  with  a 
crof  3m  and  a  time  constant  of  120  seconds.  These  values  are  based  on  the  fact  that  the 
combination  of  the  satellite  clock  and  ephemeris  errors  tend  to  exhibit  1-a  values  of  3m 
and  the  long  time  constant  reflects  the  bias-like  character  of  the  state  [31].  Modeling  the 
errors  as  random  walk  produced  unreasonable  results,  because  the  errors  never  converged 
to  a  solution.  The  initial  value  is  zero  and  the  initial  variance  is  9  meters^. 

3.3.7.  System  Measurement  Models 

This  section  describes  the  nonlinear  measurement  model,  the  linearized  partial 
derivative  matrix  H  evaluated  at  the  current  state  estimate,  and  the  covariance  of  the 
zero-mean  white  Gaussian  measurement  noise  v  vector  (the  R  matrix.)  As  a  reminder, 
the  general  form  of  the  measurement  equation  is 

z{t^ )  =  h[x(r,. ),  t,.  ]  -I-  v(t,. )  (3-24) 

where 

z(r . )  =  discrete-time  measurement  vector  at  time  ti 

h[x(t.),t;]  =  nonlinear  function  of  the  state  vector  and  time 
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v(?,) 


=  discrete  time  measurement  noise  vector  with  a  zero-mean  white 
noise  process,  which  is  independent  of  y/(t) ,  and  having  covariance 
R(t,. )  defined  by: 

£;[v(0]  =  0  (3-25) 


^R(ti)for  t,=tj 
I  0  for  t,  *  tj 


(3-26) 


3.3.7. 1.  Filter  Measurement  Model  -  GPS  Code  Measurements 


The  number  of  GPS  code  measurements  is  dependent  upon  the  number  of  satellites 
visible  to  the  receivers.  If  there  are  8  satellites,  all  visible  to  the  three  receivers  (sensor  1, 
sensor  2  and  the  bomb),  there  will  be  24  GPS  code  measurements. 

The  h  vector  models  the  incoming  measurements  z  in  terms  of  the  estimated  states. 
For  GPS  code  measurements  to  sensor  1,  the  measurement  equation  is: 


Pops  pi  ~  (i^w  ”  Vpi^  +  (z^„  —  Zpi)  +  cS  tp,  -I-  Vp^  (3-27) 

where 

^sv » ysv  ’  =  satellite  position  elements 

jCpi ,  jpi ,  Zpi  =  sensor  1  position  elements 
c<^tpj  =  receiver  clock  error 


£^[t^Pi]  =  0  and  E[Vpi(ti)Vpftj)]  =  2.6m^  secS^^ 

The  value  of  2.6m^sec  is  based  on  typical  multipath  and  noise  values  [31]. 

Recall  that  the  H  matrix  represents  the  partial  derivatives  of  the  nonlinear 
measurement  model,  evaluated  at  the  current  state  estimate 


H[t,;x(tr)] 


^[x(t),  t] 


Sx 


X  =  X(((  ) 


(3-28) 


(3-29) 


Let  e  represent  the  unit  line  of  sight  vector.  The  vector  is  described  as 
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11 

(3-30) 

(3-31) 

iysv-ypi) 

(3-32) 

-Xp,f  +{y„-ypif  +iz,,-Zpiy 

1 

(3-33) 

'sji^sv  -^Pl)  y Pl^  ^Pl) 

The  row  of  the  H  matrix  corresponding  to  the  GPS  code  measurement  at  sensor  1  is 

[«!  ^2  ^3  0  •  •  •  0  1  0  •  •  •  o] 

where  the  “1”  is  in  the  34*^  column  (corresponding  to  the  sensor  1  clock  error.) 

The  measurements  for  sensor  2  and  the  bomb  are  calculated  in  the  same  manner  as 
sensor  1.  The  difference  is  the  unit  vector  e  corresponds  to  the  appropriate  position  states 
and  the  “1”  to  clock  error  state  of  sensor  2  or  bomb. 

3.3. 7.2.  System  Measurement  Model  -  SAR  range  measurement 

For  both  the  bomb  and  the  target,  there  are  two  SAR  range  measurements-  one  from 
each  sensor.  The  measurement  equation  for  the  SAR  range  measurement  is 

=  (1  +  sf,, )  X  -x,f  +{yp,-y,f  +  Vp,  (3-34) 

where 

Xp^ ,  jpj ,  Zpj  =  sensor  1  position  coordinates 
x^,y^,z,  =  target  position  coordinates 

^/pj  =  sensor  1  scale  factor 

Vp^  =  SAR  range  noise 

jE[t»Pi]  =  0  and  E[Vp^{t)Vp^{t  +  t)]  =  8  s,&c 

This  number  is  based  from  the  work  done  in  [29,54]. 
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The  row  of  the  H  matrix  corresponding  to  the  SAR  range  measurement  is 
kl  0  ••  ^  “  ^r2  "  0  ■  ■  -  0  0  •••  O] 

where 


(3-35) 

=e^x(L  +  sfp) 

(3-36) 

'r-i  =C3X(1  +  5/p_) 

(3-37) 

11 

1 

(3-38) 

d 

(3-39) 

"  d 

II 

1 

(3-40) 

d 

d  =  —  Xj)  +  (jpi  “  jy)  +  iZpi  Zf) 

(3-41) 

The  appropriate  unit  vectors  correspond  to  the  position  states  of  the  sensor  and 
target/bomb,  and  d  corresponds  to  the  scale  factor  state. 

3.3.7.3.  System  Measurement  Model  -  SAR  Range  Rate  Measurement 

Just  as  there  are  two  SAR  range  measurements  per  sensor,  there  are  two  SAR  range 
rate  measurements.  The  measurement  equation  for  the  SAR  range  rate  is 


SAR 


P  -P 

* />1  ■*  r 


rangerate 


^|(xpl-x,y +iy  PI- +(zpi-z,y 


’(Vpi  -v,)  +  v, 


rrPl 


(3-42) 


where 


P/>i  =  yp\  Zpi  ]  =  sensor  1  position  vector 
p^  =  y,  z,]  =  target  position  vector 

Vpi=[i:pi  jpi  Zpi]  =  sensor  1  velocity  vector 
zj  =  vector  of  target  velocities 
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=  SAR  range  rate  noise 


E{v„p^  ]  =  0  and  £'[v„p,  (t,  )v^py  {t ^ )]  =  0.04^  (m^  sec) 

These  values  are  based  on  work  done  in  [23, 46]. 

The  partial  derivative  H  matrix  row  for  the  SAR  range  rate  technically  should  include 
terms  for  the  sensor  and  target/bomb  velocities  and  the  sensor  and  target/bomb  positions. 
While  the  row  (shown  below)  takes  the  velocities  into  account,  the  terms  corresponding 
to  the  position  are  so  small  that  they  are  ignored.  The  row  corresponding  to  the  SAR 
range  rate  between  sensor  1  and  the  target  is  as  follows: 

[O  •  •  •  0  e,  ^2  ^3  0  •  •  •  0  -  -  ^2  -  ^3  0  •  •  •  O] 


3.4.  Calculations  to  Find  Accuracy  in  Bomh-to-Target  Vector 

The  3-D  vector  from  the  bomb  to  the  target  is  not  explicitly  estimated  in  the  Kalman 
filter.  Nonetheless,  calculating  this  vector  is  the  primary  goal  of  the  filter.  The  estimate 
(bomb-to-target)  can  be  calculated  from  elements  in  the  state  vector.  Subtracting  the  true 
range  from  the  estimate  will  give  the  error  in  the  range  estimate.  The  range  estimates  and 
covariances  are  found  with  the  following  equations: 


Ax  =  Xg  -  .Xj.  =  [1  - 1] 


[  < 

XqXt 

T 

1  (J^  ^ 

L 

Xj  J 

-1_ 

(3-43) 


where  Ay ,  Az ,  ,  cr^  are  found  in  the  same  manner. 

The  error  in  the  estimate  is  found  with  the  following  equations: 
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(3-44) 


S  Ax  =  Ax,„  -  Ax„,  =  Ax,„  -  ) 

S  Ay  =  Ay„,  -  Ay„.  =  Ay„  -  -  y,__^ ) 

SAz  =  Az„,  -  Az„,  =  Az„,  -  -  torj 


?flr  j,  -  DTED  estimate  of  the  target  height 
Since  the  target’s  height  is  not  estimated,  the  DTED  value  is  used.  Also  because  the 
DTED  value  is  deterministic,  the  covariance  is  only  a  function  of  the  bomb’s  covariance. 
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4.  Simulation  Results 


4.1.  Overview 

This  chapter  presents  the  results  and  analysis  of  the  system  Kalman  filter 
implemented  in  MATLAB.  The  first  section  explains  how  the  Kalman  filter  was 
implemented.  The  second  section  shows  typical  results  from  the  Kalman  filter  and 
addresses  filter  tuning.  The  next  section  describes  the  different  simulation  parameters 
and  how  the  combinations  of  these  parameters  produce  the  12  different  simulation  cases. 
Lastly,  the  results  from  each  case  are  presented  and  analyzed. 

4.2.  System  Kalman  Filter  Results  Setup 

As  mentioned  in  chapter  3,  the  measurements  built  from  the  true  environment  were 
used  in  the  Kalman  filter.  The  true  environment  produces  both  true  state  values  and 
measurements  that  are  used  by  the  Kalman  filter.  It  generates  the  true  states  and 
measurements  using  a  combination  of  pre-determined  parameters  (such  as  the  number  of 
satellites  to  use),  randomly  chosen  parameters  (such  as  the  tropospheric  humidity  term  - 
see  Section  3.2.4.2),  and  white  noise  (such  as  pseudorange  measurement  noise).  The 
Kalman  filter  then  processes  the  measurements,  producing  estimates  of  the  true  states. 
Subtracting  the  true  from  the  estimate  gives  the  error  in  the  estimate.  This  concept  is 
illustrated  in  Figure  15. 
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Figure  15.  Performance  evaulation  of  the  Kalman  filter 


The  objective  of  the  performance  analysis  is  to  characterize  the  error  process 
statistically  [31].  To  generate  the  statistical  information,  a  Monte  Carlo  analysis  is 
performed.  Multiple  samples  of  the  error  processes  are  generated  by  the  simulations,  and 
the  sample  statistics  are  computed  directly  [31].  For  this  research,  twenty-five  samples  of 
the  error  process  are  run  for  each  case.  From  these  runs,  the  mean  and  1-a  values  of  the 
state  values  are  computed.  In  addition  to  viewing  the  east,  north  and  up  errors 
individually,  analyzing  the  2-D  (horizontal)  and  3-D  errors  can  be  useful.  At  each  time 
epoch,  the  Distance  Root  Mean  Square  (DRMS)  is  calculated  by 


DRMS  = 


X[(^Ax,)^-K^Ay,)^] 


(4-1) 


where 

n  =  number  of  Monte  Carlo  runs 

d  Ax, ,  S  Ay,  =  X  and  y  error  values  for  Monte  Carlo  run  i 


The  DRMS  value  gives  the  root-mean-square  of  the  horizontal  errors.  If  the  error 
statistics  in  the  x  and  y  directions  are  similar,  then  approximately  63%  of  the  horizontal 
errors  will  be  less  than  the  DRMS  value. 
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The  3-D  version  of  the  DRMS  is  the  Mean  Radial  Spherical  Error  (MRSE).  This 


error  is  calculated  by 


MRSE^ 


^^[(SAx,f+(S^y,f+(SAz,)^ 


n 


(4-2) 


where 

n  =  number  of  Monte  Carlo  runs 

S  Ax j,S  Ay^ ,S  AZj  =  x  y,  and  z  error  values  for  Monte  Carlo  run  i 


The  MRSE  value  gives  the  root-mean-square  of  the  spherical  errors.  If  the  error  statistics 
in  the  x,  y  and  z  directions  are  similar,  then  approximately  61%  of  the  total  3-D  errors 
will  be  less  than  the  MRSE  value. 


4.3.  Filter  Output  and  Tuning 

To  give  an  example  of  the  type  of  results  the  system  filter  produces,  the  next  group  of 
figures  give  the  East-North-Up  position  error  statistics  for  the  bomb  and  the  target  for  an 
example  where  all  in  view  satellites  are  used,  the  satellite  common  errors  are  estimated 
and  the  bomb  is  tracked  with  the  same  sensor  as  the  target. 

Figure  16  shows  the  error  in  the  bomb  position  for  the  east  direction.  For  the  top 
graph,  the  multiple  dashed  lines  are  the  errors  from  25  Monte  Carlo  runs.  The  solid  lines 
represent  the  filter-computed  standard  deviations  of  the  Monte  Carlo  runs.  In  the  lower 
plot,  the  dashed  line  represents  the  mean  error  of  the  Monte  Carlo  runs,  the  solid  lines 
represents  the  error  standard  deviation  of  the  Monte  Carlo  runs,  and  the  dotted  lines  are 
the  filter-computed  l-cr  error  value  from  the  runs.  The  multiple  covariance  lines  are  due 
to  the  varying  satellite  geometry  (number  of  satellites  and  locations). 
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From  the  top  graph  in  Figure  16,  the  runs  show  little  noise  in  each  run  up  to  about 
620  seconds.  At  this  point,  the  two  radar  sensors  are  tracking  the  bomb.  These 
additional  range  measurements  are  the  cause  of  the  slight  increase  of  noise  in  the  error. 
The  filter  standard  deviation  matches  very  well  with  the  Monte  Carlo  simulations.  The 
second  graph  shows  that  there  is  an  approximate  -0.7  m  bias  error  in  the  estimate.  This 
bias  may  be  due  to  a  lack  of  observability.  With  no  true  reference  point,  there  is  no  basis 
on  which  to  compare  relative  positions. 


Time  (sec) 

Figure  16.  Error  in  Bomb  Position  -East 


Figure  17  shows  the  errors  in  the  bomb  position  in  the  north  direction.  The  top  graph 
in  Figure  17  shows  a  distinct  drop  in  error  around  375  seconds  for  a  few  of  the  runs.  This 
is  caused  by  the  loss  of  a  satellite  from  view.  In  the  second  graph,  the  standard  deviation 
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of  the  errors  also  reflects  the  drop.  Again,  at  time  620  seconds,  the  growth  in  noise  is  due 
to  the  added  sensor  measurements.  A  positive  0.4  bias  appeared  in  the  mean  error  of 
graph  two. 


Time  (sec) 


Figure  17.  Error  in  Bomb  Position  -  North 

Figure  18  shows  the  plots  for  the  error  in  the  bomb  position  in  the  up  direction.  The 
data  for  the  up  direction  in  Figure  18  shows  the  same  satellite  drop  around  375  seconds. 
Compared  to  the  east  and  north  direction  graphs,  there  is  more  noise  in  the  error,  and  it  is 
consistent  the  entire  time.  Also,  the  constant  bias  that  was  evident  in  the  other  two 
directions  is  not  a  factor  here.  Again,  the  filter  covariance  matched  well  with  the  Monte 
Carlo  results. 
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Figure  18.  Error  in  Bomb  Position  -  Up 

Figure  19  shows  the  target  position  error  for  the  east  direction.  Comparing  these 
graphs  to  the  position  graphs  of  the  bomb,  there  are  several  key  differences.  First,  the 
position  of  the  target  comes  only  from  the  range  measurements  of  the  sensors,  while  the 
bomb  positions  are  based  on  GPS  (and  sometimes  range)  measurements.  The  shape  of 
the  filter-computed  covariance  is  a  result  of  the  changing  geometry  between  the  target 
and  the  sensors.  The  covariance  drop  at  620  seconds  is  due  to  the  pickup  of  the  bomb  by 
the  radar  sensors,  and  the  drop  at  720  seconds  is  from  including  the  atmospheric  error 
(scale  factor)  in  the  bomb  measurement  model,  which  improves  the  sensor-to-target 
measurement  accuracy.  The  two  large  changes  in  the  mean  are  due  to  the  turns  of  the 
sensors  platforms. 
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Figure  19.  Error  in  Target  Position  -  East 


Figure  20  shows  the  error  in  the  target  position  for  the  north  direction,  and  some 
characteristics  are  similar  to  the  north  target  position  error.  The  two  drops  in  the 
covariance  are  again  due  to  the  bomb  picking  up  range  measurements  from  the  sensors 
(620  seconds)  and  the  bomb  atmosphere  error  estimated  as  the  target  atmosphere  error 
(720  seconds).  The  steady  growth  of  the  error  between  230  seconds  and  600  seconds 
could  be  due  to  a  combination  of  the  modeling  of  the  atmosphere  error  as  a  scale  bias  and 
geometry  of  the  sensors.  As  seen  in  Figure  21,  the  true  atmosphere  scale  factor  is  neither 
constant  nor  linear.  Recommendations  to  change  the  model  are  addressed  in  the 
concluding  chapter. 
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Figure  20.  Error  in  Target  Position  -  North 
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Atmosphere  scale  factor  of  Sensor  1  to  Target 
(Atmosphere  error/Range  from  sensor  1  to  target) 


Figure  21.  Atmosphere  scale  factor  of  Sensor  1  to  Target 
(Atmosphere  error/Range  from  sensor  1  to  target) 

Note  that  the  states  for  all  the  positions  are  fairly  well  tuned  (i.e.,  filter  and  Monte 
Carlo  1-a  values  are  similar).  However  this  is  not  the  case  for  all  the  states.  Figure  22 
gives  an  example  of  the  east  acceleration  state  of  the  bomb.  The  covariance  is  much 
greater  than  the  1-a  values.  This  is  not  necessarily  bad  tuning.  In  an  operational  filter, 
the  acceleration  state  process  noise  is  often  set  conservatively  high  to  ensure  that  any 
unexpected  acceleration  can  be  tracked  by  the  filter  [31].  However,  to  get  more  optimal 
performance  of  the  filter,  all  states  should  be  correctly  tuned. 
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Figure  22.  Error  in  Target  Acceleration  -  East 


4.4.  Filter  Parameters  and  Case  Definitions 

For  these  simulations,  there  are  three  parameters  of  concern:  1)  the  number  of  GPS 
satellites  used  to  produce  the  navigation  solutions,  2)  whether  or  not  to  estimate  the 
common  satellite  errors,  and  3)  whether  or  not  to  track  the  munition  with  the  same  radar 
sensors  that  are  tracking  the  target.  Each  of  these  parameters  will  now  be  addressed. 
4.4.1.  Number  of  GPS  satellites 

Regarding  the  number  of  satellites,  there  are  three  different  modes.  The  first  mode  is 
an  all-in-view  satellite  navigation  solution.  In  this  mode  each  receiver  uses  all  visible 
satellites  (any  satellite  above  5°  is  considered  visible).  In  the  second  mode, 
measurements  from  two  satellites  are  removed  for  each  receiver.  With  the  dynamics  of 
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the  sensor  platforms  it  is  quite  possible  that  all  visible  will  not  by  utilized  by  the  receiver. 
In  the  last  mode,  only  four  satellites  are  used  to  find  the  solution.  This  mode  simulates 
most  current  military  receivers,  which  only  track  4  satellites.  Like  military  receivers, 
only  combinations  of  4  satellites  that  result  in  a  good  satellite  geometry  (i.e.,  that  have 
low  Dilute  of  Precision,  or  DOP,  values)  are  used. 

4.4.2.  Estimating  the  SV Errors 

As  described  in  Section  2.4,  each  satellite  has  errors  (ephemeris  and  clock).  These 
errors  are  common  to  receivers  that  are  in  close  proximity.  Estimating  these  errors 
should  improve  the  position  estimates.  When  included,  the  common  satellite  errors  are 
modeled  in  states  42  through  the  end  of  the  state  vector  (see  Table  4). 

4.4.3.  Tracking  the  Bomb  and  Target 

The  benefit  of  tracking  both  the  bomb  and  the  target  with  the  same  sensor  is  that  the 
errors  of  the  SAR  range  measurements  are  correlated,  so  they  can  be  observed  and 
removed  by  tracking  the  bomb  (which  has  a  GPS-based  position).  The  SAR  range 
atmospheric  error  is  determined  by  the  amount  of  atmosphere  the  signal  must  travel 
through.  As  the  bomb  approaches  the  target,  the  amount  of  atmosphere  between  the 
sensor  and  the  bomb  begins  to  match  the  atmosphere  from  the  sensor  to  the  target. 
Therefore  the  error  in  the  sensor-to-bomb  signal  and  the  error  in  the  sensor-to-target 
signal  become  more  correlated.  Since  both  the  sensor  and  the  bomb  have  relatively 
precise,  GPS-based  positions,  the  distance  between  them  is  accurately  known  as  well. 
The  filter  can  use  this  information  to  determine  the  atmospheric  error  between  the  sensor 
and  the  bomb,  and  because  of  the  correlation,  information  about  the  sensor-to-target 
atmospheric  error  is  obtained. 
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4.4.4.  Case  Descriptions 

By  varying  the  number  of  satellites,  whether  or  not  the  satellite  common  errors  are 
modeled,  and  whether  or  not  the  bomb  is  tracked  by  the  radar  signal,  12  cases  were 
created.  Table  5  shows  each  of  the  12  cases. 


Table  5.  Case  Descriptions 


Bomb 

Tracking 

SV  Error  Model 

Case  1 

Case  5 

Case  9 

No  SV  Error 

Model 

Case  2 

Case  6 

Case  10 

No  Bomb 

Tracking 

SV  Error  Model 

Case  3 

Case  7 

Case  1 1 

No  SV  Error 

Model 

Case  4 

Case  8 

Case  12 

All  SVs 

Limited  SVs 

2  fewer  SVs  per 
receiver 

Limited  SVs 
Total  of  4  SVs 
per  receiver 

Before  the  all  results  are  summarized,  data  from  case  1  will  be  presented  to  give  an 
example  of  how  the  statistics  are  calculated. 

4.4.5.  Case  1  -  All  Satellites,  Modeling  Satellite  Error,  Bomb  Tracking 

The  following  figures  present  the  results  for  case  1  in  graphical  form.  Each  figure  has 
two  graphs:  1)  the  Monte  Carlo  runs  with  the  filter-computed  1-g  bounds,  and  2)  the 
mean  error  and  standard  deviation  from  the  Monte  Carlo  runs,  along  with  the  filter- 
computed  standard  deviation. 

Figure  23  shows  the  East  error  for  the  bomb-to-target  vector.  Looking  at  the  second 
plot,  a  large  change  in  the  mean  occurs  approximately  at  230  seconds,  620  seconds,  and 
700  seconds.  The  first  time  reflects  the  first  dynamic  turn  that  the  sensors  perform.  The 
second  time,  620  seconds,  marks  the  release  of  the  bomb.  The  third  time  marks  the  time 
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that  the  SAR  atmosphere  range  error  is  applied  to  the  bomb.  This  is  significant  because 
the  SAR  atmosphere  error  is  now  correlated  between  the  target  and  bomb,  resulting  in  a 
better  estimation.  In  previous  filter  implementations,  the  SAR  atmosphere  estimate  was 
applied  to  both  the  target  and  bomb.  However,  the  errors  in  the  SAR  atmosphere  error 
were  significantly  wrong.  When  the  SAR  atmosphere  errors  were  applied  to  the  bomb 
only  after  it  was  lower  than  an  altitude  of  4000  meters,  the  errors  were  more  reasonable. 


Figure  24  shows  the  north  errors  in  the  bomb-to-target  vector.  For  the  north  error,  the 
mean  also  changes  at  the  first  two  specified  times.  In  comparing  the  amount  east  error 
with  the  north,  the  north  error  is  less.  This  is  due  to  the  geometry  of  the  sensors  relative 
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to  the  bomb  and  target.  Recall  that  both  sensors  are  always  to  the  south  of  the  bomb  and 
target  (see  Figure  9).  Therefore,  the  range  measurements  provide  more  information  about 
the  north  than  the  east  axis. 


Figure  24.  Case  l-North  Error  Bomb-to-Target 


The  up  direction  Monte  Carlo  runs,  shown  in  Figure  24,  have  less  biased  errors  and 
much  higher  standard  deviations  than  the  horizontal  axes.  This  is  due  to  the  DTED  error, 
which  is  a  randomly  chosen  constant  for  each  run.  Had  the  DTED  error  been  ignored, 
then  the  results  would  show  less  bias  (the  1-a  DTED  error  value  is  10m.) 

The  graphs  from  the  other  1 1  cases  are  available  in  the  Appendix. 
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Figure  25.  Case  1-  Up  Error  Bomb-to-Target 


4.5.  Results  Summary 

To  summarize  the  results  of  the  simulations,  the  statistics  of  the  various  cases  have 
been  put  into  tabular  form  in  order  to  facilitate  comparison.  These  statistics  include  the 
bomb-to-target  error  means  and  standard  deviations  in  the  ENU  frame,  the  error  DRMS 
in  the  East  and  North  axis  (horizontal),  and  the  Mean  Radial  Spherical  Error  (MSRE)  in 
all  three  axes.  Because  this  is  a  bomb  application,  a  “terminal  period”  is  specified 
(defined  as  the  last  50  seconds  of  the  bomb  flight).  This  period  is  chosen  because  the 
final  approach  is  crucial  to  the  navigation  solution  of  a  bomb.  In  order  to  produce  a 
single  number  for  the  whole  terminal  period,  the  mean  errors  are  averaged  over  the  time 
period,  and  RMS  values  of  the  standard  deviations,  DRMS,  and  MRSE  errors  are 
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calculated  over  the  terminal  period.  (These  RMS  values  are  essentially  “average”  values 
for  those  statistical  measures). 

The  first  section  (below)  compares  the  four  cases  with  all-in-view  satellites  solutions. 
The  next  section  compares  the  limited  satellite  cases  where  two  satellites  have  been 
pulled  from  the  solution.  Thirdly,  the  four  satellite  cases  are  compared.  Finally,  cross 
comparisons  between  the  different  satellite  cases,  the  bomb  tracking,  and  satellite  model 
estimation  are  addressed. 

4.5.1.  All-in-View  Satellite  Cases 

The  all-in-view  satellite  cases  are  cases  1  through  4.  The  results  are  found  in  Table  7. 


Table  6.  All-in-view  Satellite  Cases  (1-4) 


Case 

sv 

used 

Bomb 

Track? 

SV 

error 

mod¬ 

eled? 

East  (m) 

North  (m) 

Up(m) 

DRMS 

(m) 

MRSE 

(m) 

mean 

std 

mean 

Std 

mean 

Std 

1 

All 

Yes 

Yes 

1.30^ 

2.38 

1.90 

rti' 

11.32 

3.45 

11.65 

2 

All 

Yes 

No 

1.17 

2.26 

-0.31 

2.04 

11.38 

3.46 

11.74 

3 

All 

No 

Yes 

i  3.75 

3.56 

-0.39 

3.18 

4'  2.33 

12.47 

13.82 

4 

All 

No 

No 

3.55 

-0,97, 

3.17 

9.72 

•  5.90 

11.27 

In  examining  the  table,  the  results  provide  some  insight  into  whether  to  estimate  the 
common  satellite  errors  and  whether  to  track  the  bomb. 

Tracking  the  bomb  gives  a  definite  increase  in  accuracy.  When  tracking  the  bomb, 
the  mean  error  decreases  from  3.75  m  to  1.30  m  in  the  east  direction  (see  case  3  to  easel). 
Similar  results  are  for  the  comparison  of  case  4  to  case  2.  Tracking  also  improves  the 
standard  deviation  values.  In  tracking  the  bomb,  the  standard  deviation  went  from  3.18 
meters  to  1.90  meters  for  the  north  direction.  The  total  DRMS  (horizontal)  decrease  from 
about  6m  to  approximately  3.5m  is  about  a  2.5  meter  jump  in  accuracy. 
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Estimating  the  satellite  errors  yields  no  overall  improvement  in  accuracy  in  this  case. 
However,  it  is  interesting  that  the  mean  errors  in  the  north  direction  where  significantly 
lower  when  estimating  the  satellite  errors.  For  example,  the  mean  error  went  from  — 
0.3051m  to  -O.0478m  from  case  2  to  case  1  and  -0.97m  to  -0.39m  from  case  4  to  case  3. 
There  is  no  obvious  reason  for  this  fact,  but  it  could  be  due  to  the  geometry  of  sensors 
and  allowing  them  to  take  full  advantage  of  what  little  improvement  estimating  the 
satellite  errors  provide.  Also  the  mean  in  the  east  are  all  positive  and  biased.  This  again 
is  the  result  of  the  geometry  of  the  sensors. 

The  up  direction  does  not  give  much  information  on  the  benefits  of  estimating  the 
satellite  errors  or  tracking  the  bomb  because  the  DETD  error  dominates  the  results.  The 
true  height  for  each  case  is  zero.  If  the  DTED  error  is  close  to  zero,  then  the  amount  of 
error  caused  by  the  DTED  estimation  is  low.  To  prove  this  point,  the  actual  DTED  errors 
of  the  four  cases  were  examined.  From  this  observation,  the  reason  that  case  4  has  the 
best  results  is  because  case  4  had  lower  DTED  errors  than  the  other  cases.  Had  the  DTED 
errors  been  forced  to  be  the  same  for  each  of  test  cases,  the  results  would  probably  be  less 
random  and  could  give  some  insight  on  the  impact  of  the  parameters  in  the  up  direction. 
However,  because  the  bomb  is  expected  to  come  in  vertically,  the  impact  of  the  error  in 
the  up  direction  is  reduced.  For  this  reason,  the  results  of  the  up  direction  and  the  MRSE 
will  not  be  included  in  future  tables. 

4.5.2.  Limited  (-2)  Satellite  Cases 

Table  7  presents  the  results  of  the  first  limited  satellite  case,  where  2  satellites  are 
removed  from  each  receiver.  The  results  in  case  2  are  similar  to  the  results  in  the  full 
satellite  cases.  The  bias  in  the  east  direction  is  again  a  product  of  the  geometry  of  the 
sensor  platforms. 
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Table  7.  Limited  (-2)  Satellite  Cases  (5-8) 


Case 

SVs 

used 

Bomb 

Track? 

sv 

error 

mod¬ 

eled? 

East  (m) 

North  (m) 

DRMS 

(m) 

mean 

std 

mean 

Std 

5 

-2 

Yes 

Yes 

2.26 

''' 

1.94 

'  "  3.37' 

6 

-2 

Yes 

No 

Hi'* 

r/'-tsk 

2.62 

1.94 

7 

-2 

No 

Yes 

3.50 

3.54 

8 

-2 

No 

No 

BI1B8 

3.10 

IlSMi 

3.75 

:  6.39 

4.5.3.  Limited  4  Satellite  Cases 

Table  8  presents  the  results  of  the  cases  where  only  4  satellites  measurements  are 
used  by  each  receiver. 


Table  8.  Limited  (only  4)  Satellite  Cases  (9-12) 


Case 

SVs 

used 

Bomb 

Track? 

SV 

error 

mod¬ 

eled? 

East  (m) 

North  (m) 

DRMS 

(m) 

mean 

Std 

Mean 

Std 

9 

4 

Yes 

Yes 

3.10 

3.23 

10 

4 

Yes 

No 

3.50 

-0.71 

3.45 

11 

4 

No 

Yes 

1-76 

3.74 

3.64 

12 

4 

No 

No 

10.01 

9.36 

15.07 

While  the  first  two  satellite  sets  were  pretty  much  the  same,  the  four  satellite  set  results 
have  a  few  significant  differences.  The  improvement  in  accuracy  for  tracking  the  bomb 
is  not  as  good,  going  5.63m  to  5.08m  DRMS.  The  overall  accuracy  is  along  the  same 
lines  as  the  all-in-view  and  minus  2  satellite  cases.  The  results  of  case  12  are  not  in  line 
with  the  other  results.  All  cases  on  this  table  used  the  same  data  sets,  so  the  only 
difference  is  how  the  filter  is  implemented.  This  would  imply  that  there  is  a  significant 
benefit  in  modeling  the  satellite  errors  in  the  4-satellite  case,  especially  when  the  bomb  is 
not  tracked.  However,  looking  at  the  graphs  of  results  of  case  12  in  the  Appendix,  there 
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seems  to  be  only  two  or  three  runs  that  are  out-of-family.  Further  analysis,  to  include 
more  runs  would  be  necessary  before  a  definite  conclusion  could  be  given. 

4.5.4.  Case  Summary 

Table  10  shows  the  results  for  all  cases.  The  cases  are  organized  according  to 
whether  the  bomb  is  tracked  and  the  satellite  common  error  is  estimated. 


Table  9.  All  Cases  Summary 


Case 

SVs 

used 

Bomb 

Track 

? 

sv 

error 

mod¬ 

eled? 

East  (m) 

North  (m) 

DRMS 

(m) 

mean 

Std 

mean 

Std 

1 

All 

Yes 

Yes 

■■  ■ 

2.38 

lii'Aos" 

1.90 

'  3.45 

5 

-2 

Yes 

Yes 

V”:;.  1.24 

2.26 

lllliii 

1.94 

^4p'3.37 

9 

4 

Yes 

Yes 

j  ^07 

3,10 

i  •'0,57' 

3.23 

:Si?5.08 

2 

All 

Yes 

No 

2.26 

i  -0.31- 

2.04 

te'3.4# 

6 

-2 

Yes 

No 

— 

2.62 

llt'ilB'l 

1.94 

10 

4 

Yes 

No 

3,50 

wifli 

3.45 

,%M>'5.32 

3 

All 

No 

Yes 

^  -0.39. 

3.18 

7 

-2 

No 

Yes 

3.37 

iil^^ 

■^,,,".-0.68 

3.54 

11.’  s;98l 

11 

4 

No 

Yes 

|::.:.-b.86 

®:5.83> 

4 

All 

No 

No 

3.55 

r>0.97 

3.17 

8 

-2 

No 

No 

I- 

3.10 

llliiss 

3.75 

12 

4 

No 

No 

10.01 

9.36 

15.07' 

Overall,  the  tracking  of  the  bomb  is  the  parameter  that  most  improves  the  accuracy  of 
the  system.  As  the  bomb  is  tracked,  using  the  all  or  all  but  two  satellites  generally 
produces  the  same  level  of  accuracy.  Both  modes  are  an  improvement  to  the  four 
satellite  case;  approximately  1.2  meters  DRMS. 

The  satellite  error  estimation  typically  generally  does  not  improve  the  accuracy.  A 
vast  improvement  of  approximately  9  meters  occurred  with  four  satellites  and  no  bomb 
tracking.  To  confirm  this,  more  Monte  Carlo  runs  should  be  made  to  see  if  there  would 
be  a  change  in  the  value. 
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There  is  a  positive  bias  in  the  east  direction  mean  error  values.  This  may  be  caused 
by  a  lack  in  observability  issue  of  the  states.  In  this  system  filter,  there  is  no  “absolute 
position”  which  is  different  than  the  typical  differential  GPS  setup  with  a  “known” 
reference  system.  One  way  to  combat  this  problem  would  be  to  declare  one  of  the 
receivers  to  be  the  “reference,”  and  all  other  positions  are  based  from  its  position.  Also  a 
residual  analysis  could  be  conducted  to  determine  if  estimating  the  errors  decreases  the 
other  residuals.  As  it  stands,  estimating  the  errors  does  not  yield  much  benefit  (expect  in 
the  4  satellite  case),  and  in  a  real-time  system,  would  add  processing  time  to  the 
calculations. 

The  amount  of  error  in  the  up  direction  is  sporadic.  This  is  due  to  the  poor 
implementation  of  the  DTED  errors.  A  better  model  could  be  used.  Actual  DTED  data 
could  be  used  and  it  would  eliminate  the  entire  issue. 

Since  the  north  axis  is  more  precise  than  the  east  axis  (due  to  the  geometry  of  the 
sensors)  changing  the  location  of  one  or  both  of  the  sensors  could  increase  the  accuracy 
for  both  directions. 

In  summary,  it  has  been  shown  that  tracking  the  bomb  with  the  same  sensors  that  are 
tracking  the  target  has  a  significant  benefit  to  reducing  the  overall  DRMS  error.  When 
comparing  the  all-in-view  and  the  all  but  two  satellites  cases,  there  is  no  improvement. 
However,  the  four  satellite  cases  had  a  drop  in  accuracy  compared  to  the  first  two. 
Estimating  the  satellite  errors  showed  no  real  value  except  in  the  case  of  four  satellites 
and  not  tracking  the  bomb,  but  these  results  need  further  analysis. 
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5.  Conclusions  and  Recommendations 


5.1.  Conclusions 

Up  to  this  point  in  time,  most  of  the  error  analysis  performed  for  the  AMSTE 
program  has  been  at  the  error  variance  level,  generating  root-sum-square  (RSS)  total 
errors  from  error  budgets  consisting  of  constant  error  variances.  In  reality,  the  level  of 
error  for  GPS  positioning  and  targeting  systems  is  highly  dependent  upon  the  current 
environment  in  which  it  is  employed  (i.e.,  the  distance  between  sensor  and  target, 
altitude,  time  of  day,  etc.).  This  thesis  presented  a  comprehensive  study  of  GPS  errors 
and  their  effect  on  a  differential  GPS  bombing  system.  Satellite  clock,  ephemeris, 
receiver  clock,  troposphere,  ionosphere,  receiver  noise  and  multipath  errors  were  all 
addressed.  Ionospheric  effects  were  not  a  concern,  with  the  assumption  that  is  the 
receiver  uses  ionospheric-free  measurements.  Of  all  of  these  errors,  the  most  critical  for 
the  AMSTE  scenario  are  on  differential  tropospheric  errors  and  multipath.  Differential 
tropospheric  error  causes  problems  because,  even  after  a  model  is  applied,  the  unmodeled 
tropospheric  error  is  unpredictable  and  can  increase  as  the  baseline  and  height  difference 
between  the  receivers  increase.  Regarding  multipath,  between  20  to  60  cm  (la)  of 
pseudorange  multipath  error  is  common  for  aircraft  in  flight.  Also,  multipath  on  a 
munition  should  be  10  cm  or  less  due  to  the  limited  surface  from  whence  the  signal  may 
bounce. 

An  environment  model  was  developed  which  consisted  of  two  SAR  radar  sensors 
with  GPS  receivers,  a  falling  munition  with  a  GPS  receiver,  the  visible  GPS  satellite 
constellation,  and  a  moving  target.  The  environment  used  the  true  ranges  and  rates  of  all 
the  satellites,  sensor  platforms,  munition  and  target.  Appropriate  errors  were  added  to  the 
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true  ranges  and  rates  to  generate  simulated  measurements.  These  simulated 
measurements  will  be  used  by  the  Kalman  filter,  to  estimate  the  position  and  velocities  of 
the  sensors,  bomb,  and  target. 

Twelve  cases  were  run  (25  Monte  Carlo  runs  each)  in  the  Kalman  filter.  These  cases 
involved  variations  in  the  number  of  satellites  used  in  the  position  solution,  the  estimation 
of  common  satellite  errors,  and  the  tracking  of  the  bomb  to  the  target  with  the  same  radar 
sensor.  The  all-in-view  or  the  minus  2  satellite  mode  did  not  change  the  accuracy. 
However,  the  four  satellite  case  produced  results  that  were  not  as  precise  or  accurate. 

Overall,  estimating  the  common  satellite  errors  did  not  provide  significant  increase  in 
accuracy.  However,  in  the  case  of  4  satellites  and  no  bomb  tracking,  the  estimation  of 
the  satellite  errors  decreased  the  DRMS  accuracy  by  10  meters. 

The  most  beneficial  mode  is  the  tracking  of  the  bomb  with  the  same  sensor  as  the 
target.  In  most  cases,  the  horizontal  DRMS  accuracy  improved  from  about  6  meters  to 
approximately  3.5  meters,  just  by  tracking  the  bomb.  This  improvement  was  due  to 
modeling  the  bomb  radar  atmosphere  error  with  the  target  atmosphere  error. 

5.2.  Recommendations 

While  most  of  the  GPS  models  in  this  thesis  have  been  used  extensively  [22,36,50], 
the  SAR  model  used  in  this  research  has  only  been  used  on  three  occasions  [29,54]. 

The  following  recommendation  are  provided  to  extend  this  research: 

1.  Use  a  more  realistic  SAR  model.  The  results  of  the  simulations  in  this  thesis  should 
not  be  used  in  absolute  number  terms  because  the  SAR  model  is  not  analytically 
correct. 
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2.  Enhance  Kalman  filter  to  use  more  advanced  differential  setups.  Carrier  smoothed 
code  and  carrier-phase  measurements  are  more  accurate  forms  of  DGPS  and  could 
increase  the  accuracy. 

3.  Include  different  targeting  sensors.  Using  laser-guided  sensors  or  other  combinations 
of  sensors  could  prove  beneficial. 

4.  Change  the  flight  profiles  and  number  of  sensors  to  see  it  affects  the  accuracy.  By 
simply  moving  the  location  of  the  sensors,  a  determination  on  the  effect  of  sensor 
geometry  could  be  explored.  Additionally,  changing  the  number  of  sensors  and  the 
rate  at  which  the  sensors  provide  data  are  ways  to  improve  the  efficiency  of  the 
resources.  A  GPS  jamming  scenario  where  the  bomb  looses  GPS  and  only  sensor 
measurements  are  used  could  also  prove  beneficial. 

5.  Track  the  bomb  with  only  one  sensor  and  vary  the  tracking  rates. 

6.  Include  INS  data  with  the  GPS  measurements.  This  would  allow  for  a  large  jamming 
scenario.  Also,  different  integration  techniques  (loose,  tight,  ultra-tight)  could  be 
explored. 

7.  Increase  the  complexity  of  the  target.  For  this  research,  a  basic  model  was  used.  A 
more  advanced  model  of  the  target  could  include  stops  or  continuity  jumps  in  the 
sensor  data.  The  length  and  times  of  these  jumps  could  effect  the  accuracy. 

8.  Tune  the  Kalman  filter  to  obtain  optimal  performance.  A  single  scenario  could  be  run 
under  different  Q  and  R  values  to  find  the  best  accuracy. 

9.  Add  more  Monte  Carlo  runs.  You  can  never  have  too  many  Monte  Carlo  runs. 

10.  Use  real  DTED  data.  Real  DTED  data  would  ensure  that  the  effects  of  DTED  errors 
are  accurately  reflected  in  the  results. 
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11.  Designate  one  of  the  receivers  as  the  “master  receiver”.  The  master  receiver  would 
have  all  other  elements  referenced  to  it.  This  may  reduce  the  observability  concerns 
that  effected  the  usefulness  in  estimating  the  satellite  errors. 
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East  error  bomb-to-target  (m)  East  error  bomb-to-target  (m) 


Appendix  Graphs  of  the  Bomb-to-Target  (Cases  2-12) 

This  appendix  contains  the  graphs  of  the  other  1 1  cases  for  the  bomb-to-target  vector. 
Case  1  is  shown  in  Section  4.4.5. 
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Figure  26.  Case  2  -  East  Error  Bomb-to-Target 
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Figure  29.  Case  3-  East  Error  Bomb-to-Target 
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Figure  30.  Case  3-  North  Error  Bomb-to-Target 
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Figure  32.  Case  4-  East  Error  Bomb-to-Target 
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Figure  33.  Case  4-  North  Error  Bomb-to-Target 
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Figure  34.  Case  4-  Up  Error  Bomb-to-Target 
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Figure  35.  Case  5  -  East  Error  Bomb-to-Target 
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Figure  36.  Case  5  -  North  Error  Bomb-to-Target 
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Figure  37.  Case  5  -  Up  Error  Bomb-to-Target 


East  error  bomb-to-target  (m)  East  error  bomb-to-target  (m) 


20 
10 
0 

-10 
-20 

0  100  200  300  400  500  600  700 

Time  (sec) 

Figure  38.  Case  6  -  East  Error  Bomb-to-Target 
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Figure  39.  Case  6  -  North  Error  Bomb-to-Target 
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Figure  41.  Case  7  -  East  Error  Bomb-to-Target 
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Figure  44.  Case  8  -  East  Error  Bomb-to-Target 
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Figure  45.  Case  8  -  North  Error  Bomb-to-Target 
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Figure  46.  Case  8  -  Up  Error  Bomb-to-Target 
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Figure  49.  Case  9  -  Up  Error  Bomb-to-Target 
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Figure  50.  Case  10  -  East  Error  Bomb-to-Target 
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Figure  53.  Case  11  -  East  Error  Bomb-to-Target 
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Figure  55.  Case  11  -  Up  Error  Bomb-to-Target 
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Figure  57.  Case  12  -  North  Error  Bomb-to-Target 
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target.  The  horizontal  DRMS  position  error  during  the  terminal  phase  of  the  bomb  trajectory  drops  from  about  6 
meters  to  about  3.5  meters. _ _ _ _ _ 
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